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ABSTRACT 



An algorithm pertinent to pattern classification and 
cluster analysis of n-dimensional (n>2) data is presented. 

The algorithm simplifies the dimensionality of the sample set 
by a perpendicular projection of each point onto a plane which 
is determined by two orthonormal vectors. The operator views 
the plane on a computer graphics screen, selects a point on 
the screen with a lightpen, and attempts to achieve separation 
of classes by choosing a "direction of movement." Each itera- 
tion, a movement of a single point in a certain direction, is 
actually a reorientation of the plane into which the data 
points have been projected. These iterations continue until 
the operator decides successful separation has been accomplished. 
The method is neither restricted by statistical distribution 
assumptions nor the dimensions of the sample set. From the 
investigation of seven case studies, this algorithm emerges 
as a tremendous assistance in the analysis of multidimensional 
data . 
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I. INTRODUCTION 



The algorithm presented in this paper is an improvement 
on the separation algorithm presented in the thesis by Gilbert 
Lauzon [1]. He developed a separation algorithm to use in 
pattern classification and cluster analysis of n-dimensional 
data . 

The algorithm was written in Fortran IV for the XDS 9300 
computer installation at the Naval Postgraduate School. The 
Adage Graphics and Text (AGT) system is linked directly to 
the digital computer. The Fortran callable graphics sub- 
routines in the XDS 9300 allow for communications between the 
two systems. Also, by using the library program "MAD" of the 
AGT system, the lightpen and sixteen function switches permit 
an alteration of the display on the graphics screen. 

The computer program projects each point into a plane 
which decreases the dimensionality of the analysis from the 
original n (>2) to 2. The points on the plane are presnted 
to the operator on a cathode ray tube. The user reorients 
the plane on each iteration by selecting a point and choosing 
a "direction of movement." Using the XDS 9300/AGT-10 system, 
the operator makes the decision whether to proceed with a 
different orientation of the plane or terminate the analysis. 
The human decision is based on the pattern of points presented 
on the screen. The interaction of the man-machine system 
is emphasized. 
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Many classification problems can be approached in a 
graphical diagram. A two-dimensional plot of the features is 
drawn. From this, the human must suggest a rule that partitions 
the feature space into regions, one for each category. The 
rule can be a series of straight or curved lines depending on 
the number of categories. Since the decision is based on the 
perception of the operator, the rule is quite arbitrary. 
Therefore, the analyst must insure on a proper number of 
samples in each category. In addition, the human tries to 
determine the rule which according to his subjective judgement 
minimizes the probability of error. In essence, this means 
that the probability of a random observation from a specific 
category falling outside its respective region defined by the 
rule should be made as small as possible. The preceding 
analogy is the basic characteristics that prevail in the 
separation algorithm. The graphical diagram is the display 
screen with the projections from the computations of the digital 
computer. The rule is solely determined by the operator from 
the video presentation of the graphical diagram. 

The first modification incorporated into G. Lauzon's 
algorithm was to insure that the A and B vectors were ortho- 
normal on each iteration through the computer program. 

Basically, this means that the two vectors are perpendicular 
and each vector has a length of one. 

Secondly, the existing "X-Y" video display was changed to 
a more discernible pair of symobls consisting of the "X-0" 
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scheme. This improvement greatly enhanced the ability of the 
operator to visually identify the two groups on the graphics 
screen . 

Thirdly, supplementary video information was added to the 
display on the screen. Each iteration of the algorithm 
involves the selection of a point and a "direction of movement." 
These operations are controlled by a sequence of steps 
utilizing the sixteen function switches and the lightpen. The 
function switches have been numbered left to right and top to 
bottom for easy reference. A record of the use of these 
function switches is most beneficial to the operator. So, 
the numbers of two function switches have been displayed in the 
upper left hand corner of the screen. The one number records 
the function switch corresponding to the type of symbol 
selected and the second keeps a record of the other function 
switches. The modification was incorporated to merely assist 
the operator and did not contribute to improving the separation 
al gori thm. 
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II. DISCUSSION OF THE PREVIOUS ALGORITHM 



The previous separation algorithm was presented in the 
thesis by Lt. Gilbert Lauzon [1]. The algorithm is applicable 
to use in pattern classification and cluster analysis of 
n-dimensional data. The computer graphics system projects 
each point of a data set into a plane. The project ion is 
not necessarily orthonormal. The method reduces the dimen- 
sionality of the problem from n to 2. 

First, a description of the algorithm will be presented 
for a pattern classification problem. Assume there are two 
classes of objects, P-| and p£. These two groups are comprised 
of L-dimensi onal vectors of numbers. These numbers may be 
measurements or merely some assigned values of nature. The 
goal of the analysis is to determine the existence of a 
partition of the two groups. Once a separation has been 
achieved and a discriminant function has been determined, it 
can be used to classify new objects by the rule: 

f(Z) > 0 -> Z i s a member of the set P^ ; 

f(Z) < 0 -*■ Z i s a member of the set 

The two groups are considered separable if an appropriate 
discriminant function can be specified. Otherwise, the classes 
are inseparable and no discriminant function can be assigned. 
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The L-dimens i ona 1 vectors of the measurements on objects 
from the classes and will be denoted by X and Y, 
respectively. The algorithm commences by calculating the 
average value, DAV, of each of the L-components of the data 
set . 



DAV ( K) 



( 



1 



M 



N 



M + N 



tf) ( 2 V ( I , K) + £ X ( J , K) ) ; K = 1 ,L 



1 = 1 



J = 1 



Here M and N are the number of vectors in P-j and P^, respec- 
tively. Using the L-component vector DAV, the data points are 
centered as follows: 



X 1 ( J , K) = X ( J , K ) - DAV ( K ) ; J=1,...,M and K=1,...,L 
Y ' ( I , K ) = Y ( I , K ) - DAV ( K ) ; 1 = 1 N and K=1,...,L 



Two arbitrary non-zero vectors, A and B, in Euclidean L-space 
are chosen. A linear transformation i s computed for each point 
using the two vectors. This amounts to projecting each vector 
to a point in Euclidean 2-space with coordinates 



T A)B (Z) = (A-Z, B-Z) 

L L 

where A*Z = £ a^z.,, B*Z = £ b^z,,, Z is an arbitrary point 

K= 1 K K K=1 K K 

from the X or Y set. The coordinates are scaled so the points 
can be displayed on the AGT screen. The procedure determines 
the maximum absolute value of both coordinates as shown: 
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TMAX = max (max | A » X | , max | A • Y | ) 

J I I = 1 , . . . ,N 

J = 1 . ,M 

SMAX = max (max | B * X | , max | B • Y | ) 

J I 

Each coordinate is divided by the respective maximum value. 

T a b (Z) = (A-Z/TMAX, B-Z/SMAX) 

At this point, the video display will appear depicting the 
projections of the X and Y points. The operator selects a 
point and chooses a "direction of movement." For example, 
if the tenth X point is picked and a move .right is chosen, 
a new vector, A 1 , is computed. 



A ' ( K) = A ( K) + X ( 1 0 , K ) 

The new values of the A and 
"direction of movement" and 

Pi rection of Move 
RIGHT 
LEFT 
UP 

) 

DOWN 



K = 1 , . . . ,L 

B vectors are dependent on the 
the selected point as follows: 

New A New B 

A + Z B 

A - Z B 

v 

A B + Z 

A B - Z 



The computer returns to the linear transformation calculations 
and repeats the procedure using the new vectors. 

The goal of the procedure is a separation of X's and Y's 
into disjoint areas of the plane (screen). The operator 
proceeds using the lightpen and sixteen function switches until 
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an appropriate separation is attained or discontinues the 
analysis. If the two classes have been successfully separated, 
the discriminant function f(Z) can be specified. For instance, 
if the X's and Y’s on the screen can be separated by a straight 
line, the discriminant function would be of the form: 

f(Z) = d (A-Z/TMAX) + e (B-Z/SMAX) + f 

where d, e, and f are scalars. The function need not be 
unique. The pattern cl as s i f i cation is completed once the 
discriminant function has been determined. 

The separation algorithm for the cluster analysis problem 
is very much like that in the pattern classification. Normally, 
the set of data is not broken down into two classes. The 
entire set is presented as a set of Y's. The concept of this 
method is the determination of groupings from the sample set. 

It is normally used on data where the goal of the analysis is 
to detect clusters. The final groupings of the data points 
into the appropriate number of disjoint subsets is the decision 
of the operator. Each of the subsets represents a cluster. 



10 



III. DISCUSSION OF THE NEW ALGORITHM 



As mentioned in Section I, the separation algorithm was 
modified using orthonormal vectors, A and B. This modifica- 
tion is important for both the pattern classification and the 
cluster analysis problems. The following paragraphs will 
describe in detail’ the concept of the orthonormal projection 
al gori thm. 

First, let us assume the analysis is conducted on a pattern 
classification data set. The data has been divided into N 
vectors in and M vectors in The P^ set is read by the 

computer as the Y points and the set as the X points. The 

data points have been centered by calculating the L-dimensioned 
vector, DAV , by averaging each component of the entire data 
set . 

Next, the operator must select the A and B vectors to be 
inserted into the computer memory. If the inner product 
between two vectors is zero, then the two are perpendi cul ar . 

The length of any vector is determined by taking the square 
root of the inner product with itself. For example, the length 
of a vector A, denoted by ||A||, is given by 



Now the linear transformation is performed on each point 
as foil ows : 




L 

( A • Z , B • Z ) ; A-Z = Z 3 k z K’ B * Z 



L 




K= 1 
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This transformation is an orthonormal projection of the point 
Z (an arbitrary X or Y) into the plane defined by the ortho- 
normal vectors, A and B. The two sets of projections are 
presented to the operator by the video display on the screen. 

The scaling method described in Section II is performed 
on the data points in order for them to be displayed on the 
screen.. By using only one scale factor, the same perspective 
of the two projections will be retained. This scaling term, 
VMAX , is determined as follows: 

VMAX = max ( max ( | A*X | , | B*X | ) , max ( | A • Y | , | B • Y | ) ) 

J I 

where J = 1,...,M and I = 1,...,N. The final values of the 
coordinates from the linear transformation will be 

t a,b = ( A * Z / VMAX > b-z/vmax) 

The first iteration of the algorithm commences when the 
operator picks a point on the screen and chooses a "direction 
of movement." Actually, the computations do not "move" the 
point selected but reorients the plane upon which the entire 
data set is projected. The reorientation of plane is accom- 
plished by computing two new orthonormal vectors. Suppose 
that the user has selected an arbitrary point Z and desires 
to move it right. The first step in the computation is that 
of recomputing the A vector, 

A' (K) = A ( K) + Z ( K ) ; K = 1 ,. . . ,L 
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For the projection of Z to shift horizontally to the right, 
the new vector, B ' , should be the same as B. But this is 
impossible if A 1 and B 1 are to be perpendicular. Hence, the 
projection of Z shifts also up and down. We want this vertical 
shift to be minimal. It is best to follow the remainder of 
the concept to finding such a B vector with mathematical 
expressions . 

The appropriate symbols will be defined as: 

A and B = old orthonormal vectors 
A' and B 1 = new orthonormal vectors 
I A* B | = absolute value of the inner product 

of the vectors A and B 

The new vector A 1 must be changed so it remains constant with 
a unit length. This is done by dividing the vector (A + Z) 
by its length as such: 

A' = (A + Z ) / | | (A + Z) | | 

This is the new vector which will be used in recomputation of 
the new projections. 

The new B vector will be determined using Lagrangian 
minimization. This technique is applicable because the require- 
ment to determine the minimum value of the objective function, 

| Z * ( B ' - B) | , has the two restrictions imposed on it. Now, 
to find the B' which minimizes the difference between Z - B 1 and 
Z*B, the following problem will be solved: 
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Mini mi ze 



I Z • ( B 1 - B ) | (1) 

subject to B ' *A ' = 0 (2 ) 

||B'|| 2 = 1 (3) 

Using the Lagrangian method to find the solution, let 

F = + Z-(B'-B) - X 1 ( B ' • A ' ) - X 2 (||B'|| 2 - 1) (4) 



where X.| and X^ are the Lagrangian multipliers. Taking the 
first derivative with respect to B' and setting the result 
equal to zero: 



-jj-g-i - +_ Z - X -j A 1 - 2 X 2 B 1 - 0 (5) 

Using equations (2) and (5), the value for X-j is: 

X ] = + Z-A' (6) 



This result is true if and only if X^ has a non-zero value. 
Using equations (3) and (5), the value for X^ can be determined 
as : 

X 2 - i H Z - <*•*'> A ll- 1 - (7) 



Substituting the values of X-j and X^ into equation (5), the 
final solution to the minimization problem is: 



B' 



+ (Z - ( Z • A ' ) A 
- | | (Z - (Z-A*) A 




( 8 ) 



The final step is the determination of the sign to use for the 
B' vector. By using the inner product of Z with both sides of 
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equation (8) and the second derivative of equation (4), the 
following statements can be shown to be true: 

(1) Use + sign if sign ( Z • B ' ) f sign (Z - ( Z • A ' ) ) 

(2) Use - sign if sign ( Z • B ' ) = sign (Z - ( Z * A ' ) ) 



The same reasoning can be extended to "moving" an arbitrary 
point Z to the left. However, in selecting a point Z and 
"moving" it up, the new vector B' will be determined first. 

This time the minimization problem will be: 

Minimize | Z ♦ (A ' - A) | 

Subject to B ' • A ' = 0 

lift'll 2 * i 



The following table gives a guide to the computations 
of the new vectors using an arbitrary point Z. 



Direction of Move 
RIGHT 
LEFT 
UP 

DOWN 



New A 

(A + Z ) / | | A + Z | | 
(A - Z ( / | | A - Z | | 
A* 

A* 



New B 
B* 

B* 

(B + Z ) / | | B + Z | | 
(B - Z) / | | B - Z | | 



*Denotes the use of the minimization technique. 

These new A and B vectors are used in the recomputation of the 
linear transformation T on the data set. After the projections 
have been calculated, a new display of the plane depicting the 
two subsets will appear on the screen. The operator can make 
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another selection or terminate the analysis depending on the 
presentation of the points. 

The "X-Y" letters in the old routine were changed to the 
"X-0" symbols. This was done to make the presentation easier 
for the user to discriminate the two data sets. The similarity 
between an X and Y made it perplexing when the two different 
symbols were in close proximity. The two letters are sketched 
on the screen through a set of draw and move instructions. 

These compute r instructions are stored in a vector called IPACK. 

Video information of the function switches was also 
incorporated into the computer program. Now, when the operator 
depresses one of the sixteen function switches, the number of 
the switch is displayed in the upper left corner of the screen. 
The number in the first row indicates the function switch for 
selecting either an X or a Y point. The value of these 
commands will be an 11 or 3, respectively. The remaining 
function switches will be indicated on the second row. This 
additional video display is controlled by a series of ENCODE 
and CALL TEXTO instructions in the computer program. The role 
of the function switches is to communicate between the Fortran 
graphics package "MAD" and the 9300 computer. The display of 
the function switch numbers assist the operator in efficient 
execution of the program. 
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IV. ANALYSIS OF THE SEVEN CASES 



The operation of the separation algorithm using the 
orthonormal vectors was first examined using four sets of 
artificially generated coordinates. The first set was just 
to test the modified program. The remaining three were 
generated for the sake of testing the applicability of the 
method to pattern classification. Each set of points filled 
various geometrical shapes. This quality allowed for ensuring 
proper appraisal of the procedure using mul ti -di mensi onal data 
Finally, three cases of measurements from "real world" applica 
tions were analyzed to evaluate the performance of the new 
algorithm. The seven cases that will be discussed are: 

TESTS : 

(1) One three-dimensional pyramid 

(2) Two three-dimensional cubes 

(3) Two six-dimensional cubes 

(4) One four-dimensional shell enclosing a four-space 
sphere 

APPLICATIONS : 

(5) Iris flowers 

(6 ) Fossils 

(7) Coronary heart disease 
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Prior to exemplifying on each case, a few comments are 
appropriate. The output from the computer printer comprises 
a sixty-one by thirty-seven dimensional matrix. The printout 
is a six by six inch square. The outline of the square is 
shown by asterisks and the coordinate axes are indicated by 
four zeroes. The symbols denoting the data points on the 
printout correspond to the plot on the AGT graphics screen. 

The points will appear as an X, 0, 8, E, or J. The "X" and 
"0" symbols match the positions of the X and Y projections, 
respectively. ' An "8" symbol represents two or more Y's in a 
cell; and "E" refers to two or more X's. A cell denotes a 
matrix element of the printout array. The matrix is stored 
in the computer memory. For example, PIC (31,18) in the 
computer program would refer to the symbol in the cell occupying 
row 31 and column 18 of this matrix. Lastly, the "J" letter 
indicates a mixture of X's and Y's occupying the same cell. 

The coordinate axes have been drawn in by hand. The encircled 
points signify the X or Y selected in the next iteration and 
the arrow represents the direction of movement. 

For each case, the initial values of the A and B were 
extracted from the table in Appendix IV. This table represents 
only one group of an infinite set of possibilities. This 
particular group was chosen because the absolute values in any 
dimension are nearly equal. The author felt the use of nearly 
equal components would initially bring all measurements into 
the projection by the same relative proportion. It was thought 
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that this should give a better inti ti al presentation of the 
projections than any other choice of values. However, it 
should be pointed out that the decision on the particular set 
of values to insert is left up to the user. 

Polaroid prints of the AGT plot of the projections have 
been included with the computer pr intout. The close resemblance 
between the printout and the Polaroid pictures is quite 
apparent. Of course, the Polaroid picture is a true repre- 
sentation of the projections. Any discrepancies in the plot 
from the computer printout would be due to the restrictions in 
setting up the matrix array. Adjustments of the vertical and 
horizontal controls were necessary to center the plot on the 
screen so a Polaroid print may be distorted but the relative 
distance between points has been maintained. The elaboration 
and the conclusion of each case will now be presented. 

CASE 1 

The first case involved the coordinates of a three- 
dimensional rectangular pyramid. The points represent the 
outline of this figure. The data set consisted of twenty-one 
vectors. Each consecutive iteration will show the rotation 
of the geometrical shape. The change in the A and B vectors 
merely reoriented the perspective of the figure. Also the 
shape of the figure remains true with each iteration. This 
example was used in evaluating the successful operation of the 
orthonormal projections. 
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Figures IV-1-1 through IV-1-5 in Appendix V show a five- 
step sequence in the application of the orthonormal projection 
concept. Since this case was done only as a test of the 
algorithm, no printout of the appropriate terms is given. The 
outline of this geometrical shape was drawn in by hand in the 
figures. The Polaroid picture of the AGT display, Figure 
IV-1-5P in Appendi-x V, is included to compare with the print- 
out given in Figure IV-1-5. 

CASE 2 

Two three-dimensional cubes were generated for analysis 
using the pattern classification technique. The X and Y 
subsets comprised twenty- tup 1 es depicting the outline of the 
two geometrical figures. The Y subset has its closest vertex 
to the origin at the point (3,2,5). Each of its sides has a 
length of two units. The other cube (X subset) had its closest 
vertex at (7,6,9) and also two units on a side. The two cubes 
were made disjoint. This made it feasible to specify a dis- 
criminant function depicting the partitioned subsets. Again, 
this was another test in evaluating the procedure but involved 
two disjoint figures. 

Figures I V - 2 - 1 through I V - 2 - 5 in Appendix V show a six step 
sequence in the application of the procedure on the two cubes. 
Separation was maintained throughout the sequence. The values 
of the A and B vectors have been indicated below the figures. 
The initial figure used the values from the table in Appendix 
IV. 
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A possible discriminant function in Figure I V - 2 - 5 can 
be specified as 

f(Z) = A- (Z - DAV ) /VMAX 

where A has the values printed from Figure IV-2-5. If 
f(Z) > 0, Z is a member of the subset of Y's; if f(Z) < 0, 

Z is an entity in the X's. These results can be checked by 
the following example. The printout indicated the values 
for VMAX and DAV; 

DAV ( 1 ) = 6.000 DAV ( 2 ) = 5.000 

DAV ( 3 ) = 8.000 VMAX = 4.333 

Assume that Z has the coordinates of the point Y(l) in 
Appendix III. The calculations indicate the number for f(Z) 
is + 0.662. The computation has shown that Z i's , in fact, a 
member of the subset of Y's. Also any other point from the 
Y's will satisfy the condition that f(Z) will be positive. 

One of the two cubes in Figures I V — 2 — 4 and IV-2-5 have 
been drawn in by hand. Also a figure of the Polaroid print 
is included to compare with Figure IV-2-5. 

CASE 3 

For the third application using the new separation algo- 
rithm, the dimensional space of the cube was increased from 
three to six. Seventy-tupl es were generated for each subset, 

X and Y. The unit cubes were generated on the IBM 360 computer 
as uniform distributed random numbers. A point from the Y 
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subset would be denoted by Y ( I ) = (Y(I,1), Y(I,2), Y( I ,3 ) » 

Y ( 1 , 4 ) , Y ( 1 , 5 ) , Y ( I , 6 ) ) where I = and 1 £Y(I,K) 

£ 2, K = and for the X subset X(J) = (X(J,1), X(J,2), 

X ( J , 3 ) , X ( J ,4 ) , X ( J , 5 ) , X ( J ,6 ) ) where J = 1,...,M and 0 < 

X(J,K) £ 1, K = 1,...,6. Again, the two cubes were generated 
so the two samples were disjoint. First, if one of the itera- 
tions showed a combination of "X's" and "0's" in an area of 
the plane, it would be easier to separate the two groups on 
subsequent steps. Secondly, once separati on between the two 
groups was attained, a discriminant function could be defined. 

The printout of the graphical display is shown in Figures 
I V - 3 - 1 through I V - 3 - 6 in Appendix V. The degree of complexity 
for the two samples was increased from that of Case 2, but 
a partition of the two groups was still achieved. Since this 
case is so similar to Case 2, the discriminant function will 
not be determined. Consequently, the A and B vectors will 
also be omitted. In Figure I V - 3 - 2 , the two cubes shown by 
their projections can not be separated because of the mixture 
of "X's" and "0's" around the center of the graph. It was 
easy, however, to separate the two as shown in the next step, 

Figure IV-3-3. The similarity between the printout and the 
Polaroid snapshot is shown in Figures I V - 3 - 5 and IV-3-5P. 

CASE 4 

The previous two cases were applications of the parti- 
tioning of geometrical shapes which were separable by a hyperplane. 
A sphere and a shell in four-space were created for the purpose 
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of making the division of two groups more complex. Twenty- 
tuples were generated for the X subset and forty-tuples for 
the Y subset. A typical position of a vector in the Y set is 
Y ( I ) = ( Y ( 1 , 1 ) , Y ( I ,2 ) , Y ( 1 , 3 ) , Y ( I , 4 ) ) , where I = 1 ..... N . 

The components must satisfy the requirement 2 £ Y ( I , 1 ) + Y ( I . 2 ) 

+ Y(I,3)^ + Y ( 1 , 4 ) ^ £ 4. Similarly, an X point must fulfill 
the condition that X ( J , 1 ) ^ + X(J,2)^ + X ( J , 3 ) ^ + X ( J , 4 ) ^ £ 1. 

The two subsets are again disjoint and they form a sphere 
surrounded by a shell separated by one unit in radial distance. 
Two different trials on the sample set were conducted. The 
first run was finally terminated after twenty iterations without 
achieving successful separation. It should be noted that the 
arrangement of the two groups made it impossible for an appro- 
priate partition. Of course, this is assuming that the shell 
has the proper density of points and no voids exist in it. 
However, if the procedure is modified, separation can be 
feasible. The second run was carried out by including the 
square of each component in each entity. The dimension of the 
sample set becomes twice the original. To clarify this last 

remark, if there was an arbitrary point, Z = (Z^ , 1 ^ , 1 ^ ,Z^ ) in 
2 

the original set, it will now become Z = (Z^ ,Z^ ,Z^ ,Z^ , 

2 2 2 

Z 2 ,Z^ ,Z^ ). The second trial was conducted using this eight 
component vector on both sets. 

The first trial on case four appears in Figures I V - 4 - 1 
through IV-4-5 in Appendix V. These figures confirm that the 
separation was never attained. Only the first six steps have 
been included in the figures. The remaining steps were omitted 
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since they did not improve the pattern of the projections. 

Thus, a discriminant function can not be specified. Conse- 
quently, the printout of the A and B vectors will be eliminated 
from the discussion. 

The situation improves in the second run using the eight 
component tuples. This is depicted in Figures I V - 4 - 6 through 
IV-4-10 in Appendix V. Separation was attained in a mere three 
steps of the algorithm. This modification of the procedure 
can be attempted on samples that can not be successfully 
partitioned using the basic method. The one subset may be 
enclosed inside a large sphere -like shape including the other 
subset. This is one approach to obtaining a division of the 
two groups. The components of the A and B vectors in the 
second run are printed on the bottom of each of the figures. 

CASE 5 

This sample set represents the first of three "real world" 
cases. The data set for this case portrays multiple measure- 
ments on three species of Iris flowers. The data was originally 
used by Ronald A. Fisher in his classic paper on discriminant 
analysis [3], It has become widely used in illustrating cluster 
analysis procedures. The four measurements collected are 
(1) sepal length, (2) sepal width, (3) petal length, and (4) 
petal width. The sample set has been reduced from the original 
150 to 75 samples for this particular analysis. 
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This particular sample set has the characteristics for 
either a cluster analysis or pattern classification problem. 

The author will approach the investigation using the pattern 
classification technique. It is quite possible from looking 
at the values of the three species that the one. Iris setosa, 
will have projections in an area of the plane by itself. The 
projections of the other two species appear to be nearly the 
same. Hence, it may be somewhat difficult to partition these 
latter two species. Therefore, the observations from Iris 
versicolor and Iris virginica will be divided into two different 
subsets. The Y subset will be composed of Iris setosa and 
Iris versicolor; while the X subset of Iris virginica. Each 
species of flowers was composed of twenty-five vectors in 
four dimensional space. 

The complexities in measurements is apparent in the parti- 
tioning of these three groups. Outliers may exist which 
indicate an observation is far removed from the main body of 
the data. They result from the fact that the sample under 
investigation may contain a certain proportion of extreme 
cases or from a population other than the sample from which 
the rest of the data came. An observation from the Iris 
versicolor species presents itself as a possible outlier in 
the analysis of the graphic plot. This point is Y(28) in 
Appendix III which coincides with the projection of X(24). 

The two appear as one point in the last two figures as a "J" 
symobl . The point Y ( 2 8 ) could not be separated from the group 
of X's as it is projected to the left of the line denoting 
the discriminant function. 
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Figures IV-5-1 through I V - 5 -6 in Appendix V show the 
twelve step sequence in the analysis of the three Iris 
flowers. Even though a discriminant function is generally 
determined in partitioning two groups, two linear functions 
will be specified which divide the three types of flowers. 
This will be done using Figure I V - 5 - 6 . The first functi on 
will separate the- X subset from the Y subset. Then the 
second will divide the two Iris species comprising the Y 
subset. First, the values of the A and B vectors from the 
bottom of the figure must be recorded. 

Now, a linear discriminant function has the form given in 
Section II as ; 

fm _ d ( A ( Z - DAV)) + e ( B ( Z - DAV)) + f 
TU ' " V MAX VMAX 



where d, e, and f are scalars. The f scalar is the negative 
sign of the ordinate intercept while the term d has the 
negative sign of the slope of the line. Finally, the value 
of the coefficient e will be one. Using the line in Figure 
I V - 5 - 6 separating the X and Y groups, an appropriate dis- 
criminant function has the form: 



f(Z) 



6(A(Z - DAV)) B ( Z - DAV) 2 25 
VMAX VMAX 4,40 



The magnitude for DAV and VMAX from this figure are: 
DAV(l) = 5.872 DAV ( 2 ) = 3.061 

DA V ( 3 ) = 3.804 DAV(4) = 1.212 

VMAX = 3.033 
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Using these figures, the value for f(Z) is - 0.054 when the 
components of X ( 2 4 ) from Appendix III are substituted into Z. 
Also carrying out the same computations, the result for f(Z) 
is + 0.252 for the point Y(26). These calculations have 
merely emphasized that a positive value for f(Z) implies Z 
is a member of the Y subset. If f(Z) is negative, then Z is 
an entity of the X subset. 

Now, a description of the line splitting the Y subset 
must be resolved. The form of the discriminant function has 
been determined as: 



f ( U) 



A • ( U - DAV) 
VMAX 



0.370 



where the number 0.370 is the value of the intercept on the 
abscissa axis. If the coordinates of Y(38) are substituted 
for the arbitrary point U, the result of f(U) is - 0.030. 
This calculation indicates that any point from the set of 
Y's, where f(U) < 0, represents an entity from the Iris 
versicolor species. Any U where f(U) > 0 corresponds to an 
observation from the Iris setosa. 



CASE 6 

The sample set for this analysis problem is quantitative 
measurements of morphologic parameters of specific fossils. 
The data set was extracted from Wright and Switzer's paper 
titled, "Numerical Classification Applied to Certain Jamaican 
Eocene Nummulitids" [4]. Eighty-two eight component vectors 
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were used in the investigation. Six of the original observa- 
tions had missing parameters and were therefore omitted from 
the sample set. 

Again like Case 5, this sample will be divided into two 
groups for a pattern classification problem. The author feels 
the separation will be easier to attain using two subsets. 

It should be noted that Lauzon investigated these specimens 
using cluster analysis procedures [1]. The Y subset was 
initially composed of the first thirty-eight observations. 

It is comprised of specimens from the Opercul i noi des wilcoxi 
species. This species will be denoted by Group A. The X 
subset was composed of the other three species: Opercul i noi des 

jennyi , Opercul i noi des cushmani , and Nummulites striatoreti- 
culatus. These three types of Nummulitids will be called 
Groups B, C, and D, respectively. Using these subsets, it 
was extremely difficult to pick out the projections from 
Group B and those from Group C. In order to eliminate this 
confusion, the two subsets were rearranged. Observations 
from Groups A and B were combined to form the Y subset; 

Groups C and D joined to make up the X subset. The second 
run was conducted using the above groupings. Also to make 
the investigation similar to the discussion in Wright and 
Switzer's paper, only the first six parameters were used. 

The results from the analysis were quite remarkable. Even 
on the first step, three different groups were apparent; 

Group A, Groups B and C, and Group D. In Wright and Switzer's 
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investigation, ninety-six iterations were conducted to deter- 
mine the coefficients of a linear discriminant function. They 
were able to find two different groupings. Groups A and D and 
Groups B and C. The same procedure using 96 iterations was 
repeated to find another linear function dividing Group A and 
D. At last, a linear function splitting Group B and C was 
determined. Then-the two researchers utilized these two final 
divisions as input to a standard discriminant analysis program 
The output resulted in two canonical variables. The plot of 
these two variables revealed three groupings: Group A, Groups 
B and C, and Group D. The two concluded that the specimens 
from Groups B and C could not be partitioned. The results 
in this thesis are in agreement with these findings. The 
procedure used by this author appears to be superior in two 
ways. First, there are substantially fewer computations, 
hence the computer time to accomplish the analysis is consider 
ably less. Secondly, the plot of the results is readily 
available from the printout of the computer. It should also 
be noted that G. Lauzon's investigation of this sample set 
showed only two groupings. 

Figures IV-6-1 through I V - 6 - 3 in Appendix V show the 
graphical display obtained on the first run. These figures 
depict the eight-dimensional analysis of the fossil measure- 
ments. In Figure I V - 6 - 2 , the three groups have been labeled. 
There is no visual method to distinguish between the B and C 
groups. The displays of the second run are given in Figures 



29 



IV-6-4 through I V - 6 - 7 . The values of the A and B components 
of these figures are given below each. 

The two discriminant functions separating the three 
groupings will be specified using Figure IV-6-7. First, the 
tabulation of the numbers for VMAX and DAV are as follows: 



DAV(l) = 163.220 DA V ( 2 ) = 41.43 

DAV ( 3 ) = 9.390 DAV(4) = 25.012 

DAV( 5 ) = 88.866 DAV(6) = 677.817 

VMAX . = 587.954 

A description of the line splitting Groups A and D must be 
determined. The form of the discriminant function has been 
resolved as 



f(Z) 



B - (Z - DAV) 
VMAX 



0.635 



If f(Z) > 0, then the arbitrary point Z is a member of Group 
A, B, or C. On the other hand, if f(Z) < 0, then Z is an 
entity of Group D. The equation of the discriminant function 
describing the upper line is given as: 

Win - 0 . 6 ( A • ( U - DAV)) B • ( U - DAV) 

nu ' ‘ VMAX VMAX 

U is an observation from Group A or D if f(U) > 0. Likewise, 
U is a vector from Group B or C if the result of f(U) < 0. 
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CASE 7 



The final analysis was conducted on a thirteen-dimensional 
pattern classification problem. Coronary patients from 
Letterman Hospital comprised the sixty-tuples making up the 
Y subset [1]. The other subset was composed of fifty Army 
personnel at Silas B. Hays Hospital who had undergone annual 
physical examinations. Those patients from Letterman had 
been diagnosed as having some type of coronary heart disease. 
The X subset consisted of thirteen features from a random 
sample of the Army personnel. There were no indications from 
the latter group to either classify or not as having heart 
disease. If successful separation is obtained and a discrimi- 
nant function is computed, then it could be used to classify 
new patients. 

A few of the components in each thi rteen- tupl e represented 
measurable entities such as weight and age. Other character- 
istics such as race and smoking behavior were assigned 
arbitrary values rather than numbers to signify a measurement. 
The analysis will show that these classifiers are not really 
appropriate to use. The actual measurements exhibit at least 
the characteristics of an interval scale. This is a minimum 
requirement in computing a meaningful average value of a 
parameter. Those components which had an arbitrary value 
assigned do not meet this criteria for an interval scale. In 
fact, they have the characteristics of an ordinal scale. 
Therefore, when the average values of the components, DAY, 
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are computed in the program, it is inappropriate to use these 
classifiers in the analysis. 

The first run was conducted on the thirteen-dimensional 
sample set. The displays of the projections are given in 
Figures IV-7-1 through IV-7-6. The last figure had the best 
presentation with respect to partitioning the two groups. A 
line dividing the"X's" and "0's" was drawn in the figure. 

If a discriminant function had been specified for this line, 
thirty-six of the fifty X points and forty-nine of the sixty 
Y points would have been properly classified. This result 
means that approximately 25% of the total sample have been 
projected out of their respective regions. This last state- 
ment does not seem satisfactory for any prediction model or 
se pa rati on . 

The values of the A and B vectors used by G. Lauzon of 
this data indicate a definite pattern. For those assigned 
values, the number in both vectors for each component are less 
than or equal to the absolute value of two. On the other 
hand, one of these two numbers is very much greater than 2 
for every measured component. For this reason, the components 
of assigned values will be deleted in the second run. In 
addition, this algorithm is more appropriate for separation 
of data where the sample set is on an interval scale. 

In the second run, an arbitrary observation Z will be 
denoted by 

Z ( H ) = ( Z ( H , 2 ) , Z ( H , 3 ) , Z ( H ,4 ) , Z(H ,1 2) , Z(H ,1 3) ) 
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where H = 1,...,M for an X and H = 1,...,N for a Y. Figures 
I V - 7 - 7 through IV-7-12 show the printout of the five-dimensional 
plot. The figure with the best partitioning is IV-7-12. This 
time the discriminant function would result in a smaller 
probability of error. There are fifty-four Y projections in 
the left region and thirty-six X's to the right. Now, the 
result is less than 80% of the projections fell into their 
respective areas. The author feels that this second run has 
also led to an unsatisfactory separation of the two groups. 

The reason for this conclusion is the high density of points 
near the dividing line. In addition, the distance is rather 
insignificant for the different points of the two subsets 
across this dividing line. Therefore, the author feels that 
an attempt to specify a discriminant function for this data 
would be futile. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



It can be concluded that the improvements on the algorithm 
have resulted in achieving greater success than the previous 
application. First of all, the revised algorithm provided a 
partitioning of the generated sample sets in the first four 
cases. The new separation algorithm was also applied to three 
"real world" applications. The first investigation obtained 
a satisfactory partitioning of the three Iris species. The 
three groups were divided by specifying two discriminant 
functions. The second sample set involved a more complex 
analysis problem. The observations consisted of an eight- 
dimensional array. In addition, it had been hypothesized 
that there existed four species in the data. The display for 
this analysis indicated the presence of only three species. 

The same conclusion was also reached by Wright and Switzer [4] 
in their approach to this discriminant analysis problem. 
However, their computations in reaching this result were much 
more involved than the algorithm in this paper. The investi- 
gation by G. Lauzon showed the presence of only two such 
separate groupings. From this, the improved algorithm using 
the orthonormal vectors appears to achieve a more definite 
partitioning. Finally, the employment of this method on the 
thirteen-dimensional medical data resulted in unsatisfactory 
separation. The rejection of a partition may be due to the 
use of inappropriate parameters in analyzing coronary disease. 
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All separations are really approximations. The final 
classifications must be subjected to a variety of sample sets 
for a confirmation of the discriminant function. Even after 
this has been accomplished, the presence of outliers will 
result in some finite value for the probability of error. The 
description by the discriminant function is only the optimal 
division of the groupings. 

The use of the values from the table in Appendix IV should 
be mentioned. All investigations in this paper were conducted 
choosing the initial orthonormal vectors as listed in the table 
(Appendix IV). However, if different values had been sub- 
stituted in the components, the final separation may be better. 
In the future, the operator should utilize a variety of unit 
orthonormal vectors in the attempt of achieving a successful 
parti ti on . 

This algorithm has provided a relatively easy and simple 
approach to the classification problem. It is possible that 
the merit of this method could be incorporated into a portable 
unit. A computer module designed for the appropriate computa- 
tions would provide a means for obtaining the required 
projections. A television cathode ray tube could act as the 
display system. An interface system could be developed to 
handle the input/output operations. The outstanding feature 
using a color t.v. set is the ability to display a variety 
of colors. The various subsets of the sample set could be 
linked to a respective hue. The combination of these two 
systems would make up an ideal graphi cs /computer system in 
which the algorithm would be a major feature. 
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The de ve 1 opme nt of this algorithm could be applied in an 
unlimited variety of fields. The analysis of medical and 
taxonomy samples has already been mentioned. The method is 
pertinent to any field in which there is mutivariate data to 
be separated. For instance, measurements of meteorological 
parameters may be subjected to this method of analysis. The 
motivation of the investigation could be to determine the 
existence of conditions pertaining to the development of a 
tropical storm and those present during the infancy of a 
hurricane. This algorithm may be able to- distinguish the two 
different groupings. If this can be accomplished, then the 
model can be used as a basis in predicting future storm 
a c t i v i ty . 

The final comments are directed to future changes in the 
computer routine. First, the A and B vectors should be 
punched on computer cards and read in with the data deck. 

This would eliminate the possibility of typing in an incorrect 
component. Secondly, a video presentation of the number of 
the selected point would be most helpful to the operator. 

In addition, the printout of the plot should indicate the 
selected point by a different set of letters than the present 
"X-0" scheme. These two modifications would make it possible 
to repeat the same iterations on subsequent runs. Lastly, 
it would be beneficial to the operator if a series of sample 
sets could be read. This could be arranged by adding the use 
of another function switch. The operator would now be able to 
conduct a series of different investigations without leaving 
the graphics console. 
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CENTRALIZING data P9 I NTS 
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INPUT C6MP0NENT VALUES FOR THE A AND B VECTORS 
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XXl( J)«XX1( J)+A(K)*X( J*K) 
81 XX2( J) = XX2( J)+B(K)*X( J,K) 
SUMX1*SUMX1+XX1( J) 
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D9 101 I «2# KN/ 5 
JL 3 JL+1 

FA=YY1 ( jl)/vmax 
FB = YY2(J1_)/VMAX 
IF (NPFLG*EG.O) 
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I ^ A G ( I ) = I PACK ( TA/ T3/ 0 ) 

I M AG ( I + 1 ) = I PACK ( TC/ TE/ 1 ) 
I“AG( I+2)=IPACK(TC/TF#0) 
IMAG(I+3)«IPACK(TD/TE/1) 
IMAG< I+4)*IPACK(TD/TF,C) 
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IF(LLS( ISW/7) .LT-0) G0 T0 193 
I F (LLS ( ISW/8) »LT*0) G0 T0 204 
G0 T0 117 

CHe^SlNG the direction 0F move with THE X POINT picked 





— 


















00 




— 






— 




















00 






















00 


























X 
















■ — 


CO 




X 


00 




3 
















in 




3 
















a: 


z 




3 


00 




cn 
















3 




cn 
















UJ 


3 


• — • 


CD 


t-4 




CD 
















w 




o 
















*— « 


— 


CL 


CD 






» — « 




















— 




















UI 


— 


CD 




— 
















CD 






h- 


a\ 


tH 


CO 








00 


x 


•— « 




X 




















X 




a: 




OJ 




in 








■s 


UI 


■s 


CL 




CL 


or 




















bj 


OJ 


OJ 


OJ 


OJ 










•— i 


i 


LU 


cn 


3 


3 
















CD 




















•V 






• — « 


CD 


►-« 


















CD 






CD 


CD 


CD 


<n 












V 




























. — 


X- 


X 


X 


j — 












X 


• — • 




« — i 


















i*-v 




5 < 
















t-4 


O 




w~4 


o 


•s 


tH 




















ON(D 


CD 


<n 


CD 








•* 




n 


o 


• 


x: 


o 
















3 . 


a: 


rH rH 


CO 


CD 


CD 


CD 






x 


00 


• ' 


< 




c 


< 


















■ — ■ 


U) 


• <«-« 














X 




X 


CL 




3 


CL 


















n 


•h 
















UJ 


• — t 


3 


U. 


X 


• 


3 


h~ 
















< 


% 


D O O 


o 


o 


O 






1“ 


•— « 


CL 


•— 


3 




3 


3 
















X 


3 


CL ft- 


• 


m. 


• 


• 


X 




*— « 




X 


•s 


CL 


00 


*K 


CL 
















3 


CO 


X 


\— 


X 


1 — 


X 


Z 






> 


3 


> 


X 


•s 


> 


y- 


















•—* 


3 CD 


3 


—1 


— I 


3 


•— • 




» — « 


UI 


CD 


UI 


3 


■ — ■ 


3 


3 
















• 


•s 


CD CD 


• 


• 


• 


• 


CD 




» — « 


o 




o 


(D 




O 


CD 








* 






M 


3 


> 




. — . 


— 




. — * 


CL 






»-f 




•— » 




• — * 


« 
















»-« 


Z 


UJ 


- — . ^ 


in 


vO 




00 






00 




O 




<— *■ 


or 


w 


















• 


a 


o o 


•* 


«v 


■s 


•s 


X 




in 




• 


CL 


O 


*—< 




O 








> 






X 


— * 




• • 


3: 




3 








■s 


CD 


UJ 


X 


• 


o 


X 


• 








*-A 






* — « 


3 


■w 


3 UJ 


CO 


CO 


CO 


co r-v 


< 




<r x 


zz 


a 


3 


CD 


X 


3 








D 






ti 


'*-* 


CO 


z cd 


•— « 


»— • 


*— * 


*— * T-l 








X 


• 


< 


X 


» — i 


< 










*-» 


00 


tH 


3 


JZ 3 


z 


• • 




— 


■ — 


" — *H 


CD 


CO 




UJ 


CL 


CL 


• 


- — 


or 


• 










« 


+ 




<C 3 


3 


x x 


CO 


CO 


CO 


CO 


Z 


H 


UI 


X 


UJ 


CD 


or 


O 


CD 


or 


O 








» — « 




W 


X z 




UJ CO 


3 


3 


3 


-J CD 


* — « 


X 


o 




•— ft 




3 


CD 




3 


it 


o 


—i 00 


♦ 


♦ 


3 


3 ■ 1 


3 




_J 


3 


— J 


3 X 


X 


X 


CD 


3 




3 


* — « 


X 


3 


►H 


cn 


II 


II 


in in 


fH 


3 X ll 


3 


— 


— 


— 


— ' 


>— 


O 


UJ 


CJ 


3 




3 


■ — 


— 


3 





o *- 


1— 




n 


II 




— Z Q 


< 


u_ 3 


3 


3 


3 


3 <D 


*— • 


X 


Z 


< 


3 


C 


3 


u 


*< 


3 


LU W> 


CL 


CD 


X 


rvl CD 


3X3 


u 




•— t 


— • 


•— 


•— CD 


CL 


•— 


UJ u 




LJ 


►— * 


» — < 


U 


• — * 






X o 


•— 




o 


•—3 3 


in 


























00 




















*-♦ in 


OJ 














OJ 












00 




















^ 00 ^ 



H T-* t- 4 fH « 



U 



o o u 



u u 






43 



WRITE (6,147) JED 
G9 T9 117 



ro 








— 


4 


— 


— 




LO 






X: 






— CO 


m ro 


_J 




— </) 


a: 2 : 


_j in 


CD 




oc 2 


LU U_ 


CO ^ — * 


O 


in 


UJ U_ 




CD 


3 


o 


►-< — 




- (D 


— 







co or lj ^ x ^ 
uj — cc ~ or 

*— • *v tiJ CD Cf Lu 
^ OJ *-« CD UJ •-• 



■s 3 ~ 
»-< O ~) *H 
^ »H <* O 

co — o *-« 
^ x <t — 



CVJ *—* 

«* o 

CD —i 
c — 



CD 



CD 

CD 



X 

CD 

or 

x 



ro cr 

^ UJ 



c\j o 

^ «h N 

CO — *-* 





UJ 


•— * 


3 




X 


• 


2 


X 
















• — 










CD 




LiJ 


►— € 


3 






X 


» — * 


CL 


•— « 


ID 




3 


3 
















CD 










2 




X 


►H 


CL 


CD 




»-» 


•s 


X 




CL 


00 


% 


CL 
















< 






Q 




• — i 








X 


X 




• — - 


> 


ID 


> 


X 




> 


X 
















2 






UJ 




> 






> 


3 






•— * 


UJ 


O 


Ul 


3 


— ■ 


UJ 


3 
















» — « 






3 




<D 






UJ 


(D 


CD 






Q 




o 


CD 


OJ 


o 


CD 
















• 










2 




• — i 


o 




CD 






•— « 


- — ■ 


• — « 




'■ — 


*— 1 
















3 


UJ 
















►H 








00 


— 


o 






QC 












2 






«y 


2 










X 




00 


— 


O 






in 




• 


CL 


O 


•— r 


►— < 


O 














X 


• 






v0 




00 




Ln 




• 


O 






<D 


Ul 


X 


• 


o 


X 


• 








•H 






3 


— ^ 






w 










CD 


LJ 


• 




-a- 


X 


2 


CL 


UJ 


CD 


CL 


UJ 








U 






n 


►Hi 








r- 


CD 




-* 


X 


2 


CD 


ro 




X 


• 


c 


2 


» ■ 4 


< 


2 








3 






• — « 




UJ 




vO 


t — i 


2 




— 


X 


• 


LJ 






Ul 


X 


or 


• 


■ — • 


LL 


• 










i 


+ 




CD 


3 




' — • 


t — < 


> — « 


in 




LJ 


cr 


m 


U 


UJ 


X 


UJ 


o 


or 


o 


CD 


or 


o 








3 


3 


T-l 


< 


2 








X 


n 


UJ 


X 


3 


o 


X 


o 




•— I 




UJ 


CD 




Ul 


II 


* — < 


o 


LO 


* 


* 


v0 


2 


•— * 


3 


UJ 


CD 


< 


X 


O 




►H 


Ul 


X 


<D 


3 


— 


_J 


* — i 


2 


3 




CO 


II 


11 




vO 


sC 




3 


X 


it 


X 


X 


or 


X 


CD 


3 


— 


3 


UJ 


CD 


3 




—1 


— 




_J 


— 


o 


X 


X 




ii 


tl 




— 


2 


O 






3 


UJ 


CJ 


3 






X 


2 


< 


u 


< 


u 


U_ 


<1 


u 


Ul 


CO 


or 


CD 


X 


NJ 


3 


U 


CD 


UJ 


a: 


CD 


X 


X 


2 


< 


L 


U 




UJ 


CJ 




o 


•—» 


►— « 


CJ 


►— « 


V 






o 


3 


3 


O 


*— < 


CJ 


3 


3 


o 


►H 


►— « 


UJ 


CJ 


►H 


►H 



a\ 



ro 

LO 



N LO 

vC in vO 



Ch 

vO 



cj 



o cj cj 



44 



WR I TE ( 6/ 170 ) JED 
D9 171 K'liL 
171 A(K)=A(K)+Y(JEDjK) 

ef=o.c 



173 A(K)=A(K)/SQRT(FF) 
KS = 0 

175 E = 0 • 0 
EE = 0»0 



o 

• o 



I «”4 

■ n 

liJUi 



o o 

• • 

o o 

• • 

• • 

UiUi 



o o 
z z 
< < 



o o 

• • 

o o 



—I o 

• • 

o o 



or or 

CD CD 



o o 

«. • 

o o 



vO 



or 



ro 























• • 












U 




rH 






















ft- ft- 












u 


























CD CD 












_J 




OJ 






















• c 






r- 
































Ld Ld 






00 






>- 




x m 






















"O— ► 
















x •* 






















• • 




u 








CD 




Ld — 


V 






C 














O O 






CD 






z 




ft— *— • 


- 






♦ 














z z 




ft- 


ft- 






., * — * 






C V 




u 














< < 




or 








> 




— > 


* 






u 










CVJ 




• • 




<3 


(D 






<D 




•— U 




CO 




ft 










* 




. . 




CO 


CD 










— o 


_J V 


♦ 














* 




o o 


_i 


\ 












■% ►— « 






















• • 




- — . 








>- 




00 — 


— » o 


V 


r-» 








*4 CVJ 






0 

0 




V -4 






CD 




in 


a Ld 




U 


o 






u 


* 


** 




• * 


n 


• — 


• 










^ CD 


V O V 


Ld 






V 


* 


o cu 


hhYUO 




LO 


CD 




ft- 




LJ 




~> 










Ld 


* 


CD _J 




♦ 


Id 




vO 


Z 




— X 


>- 






— 










“> 


♦ 


• » 




Ld 


• 




CO 




o 


Ld 


NO + 


■*— 




>- 






00 


w- 


— 


U CD CD Ox 


It 


CO 


c 




I— 


u 


u *~ 


Ld 


>- 




tl 


0 

0 

7 


CJ 


X 


Ld 


w — 


rv 




V 


ti 


CD 


< 


ft- 


Q 


«—• Ld 


+ 


r~4 


— - 


• 


• 




+ 


+ 


1 


— - ^ 


tH 


* 




CO ft- 


or 


X 


CD _J 


0 


Ld 






o 


O 




u 


o 


CD 






") 




o 




u 


Ld U J 


<D Ld 


u 


<D 


— ' 


n 


it 


<D 


11 


u 


u 


u u 


<n 


— - 


u 


Ld 


CD 


x 


H- Z < 


O Ui 


Ui 


QUL 


CD O 


u 


CD 


CD 


* — » 


o CO 


*—» 


V 


CD 


* — i 


* — « 


U CJ 




vO 














00 








CTs 










T-t 










i\ 














* 














00 






<*-4 




t-4 










T~« 








*— 4 










<*-4 





uuu 



45 



IF (IER.NE.O) SUTPUT(lOl) lERdFNS 
I F ( JED • EQ • 0 ) G3 T8 117 
WRITE(6#182) JED 
D9 183 K = 1 / L 
183 A(K)=A(K)-Y(JED/K> 



— cn 
QT z 
LlJ Ll 



% «s 

ro or 

Id 






















0l 






«s 


«— i 








































z> 






OJ 


o 
















































*—4 








































V 




h- 


CO 


— ' 


«-• 










































X 


«% 


»— 


tH 






































o 




Id 




z> 








































z: 




1- 


►— » 


CL 


CD 






— - 






















u. 










•— * 




• — i 




♦“* 


l — 


O 




y 








Ll 














Ll 










> 






> 


3 




Ld 












Ll 














■ — 










CD 




•— « 


Ld 


CD 


CD 


3 




o 








' — 














h- 










51 




* 


O 




O 






Ld 








K- 








_I 


(\J 


_j or 
















*— * 








3 


3 




OJ 


3 


or 






3 




♦ 


•s 


3 










>- 




00 




O 




■ — - 








* 




3 








w— 1 


* 


— » CO 










CD 




in 




• 


O 






> 




* 


H 


CO 






*-* 


u 


— 


II 


X 
















CD 


Ld 


• 


CD 


it 


+ 


11 




II 


X 






n 


y 


y 


y 






in 




in 


CD 






h~ 




o 




y 




y 


y 


y 








y 




■ — • 




y 








vO 








X 


* 


Ld 


•k 




y 




*— > 




y 










< 












OJ 








Id 


QC 


• 


vO 








CO 












o ^r 


+ 


LD 


< 






OJ 




b~ 


n 


id 


H* 


Ld 


O 




in 


CfiOvO 


+ 




CD 






O CD 


• 00 


u 


00 


ii 




CD 


ll 


CD 


C 




o 




*— • 


Ld 


Ld 


0\ 


11 


• CD 


Ll 


CD 


ll 




o 


• CD 


O *-4 


u. 


«-• 






h- CO 


i— 


cr 


X CD 


_J 


— 


3 


b- 


r-l 


— 


O <rH 


Ll 


rH 




o 


• 


O <r“4 


II 


n 




y 


11 




o 




Ld 


Ld 


o 












y 


II 


ti 




y 


11 


o 


It 


Ll CD 


U. 


CD 




CO CD 


Ld CD 


1— 


f- 




< 


Ll 


Ll 


CL 


CD 


— • 


U CD 


Ll 


cn 


■ — 


CO 


u 


Ld CD 


b. O 


Ll. 


O 


c 


y 


CD 


y 


o 




*— « 


Ld 


U 


►— * 




3: 


O CD 


Ll O 


Ll 


o 


CO 


y 


Lu LlJ O 








in 












ro 














in 




vO 








00 






00 




00 






00 






CD 














CD 




CD 




CD 




CD 
















« — i 
























<r-4 




fH 




r-4 





y 

CO 

* 

y 

o 

LlJ 

3 

>- 

+ 

Ld 

Ld 

u 

Ld 

Ld 



cj o o 



46 



199 E=E+Y ( JED* K ) *A ( K ) 

D© 200 K= 1 / 1_ 

200 C(K)=Y(sJED/K)-EE*B(K) 



o 

• o 



I 

■ « 

Ul Ld 



o o 

• • 

o o 

• • 

t— K- 
-I 3 
• • 
Ul Ui 



o o 

7. Z 

c c 



o o 

• • 

o o 



3 o 

• • 

cd o 



a: a: 

CD CD 



O O 



00 
— CO 

a: z 

Ui u 
•— « » 
% 

r> a: 

ui 











o o 








z 




«s 




























• «-■ 








3 
































1 — 1 — 








(D 






tH 


























CD CD 








O 




OJ 


O 


























• • 


NO 






































Ld Ul 


tH 






X 




X 00 




tH 


























OJ 










X * 


X 


^ — % 
























• » 








CD 




Ul — 


3 


























o o 


U <D 






Z 




X — 


CL 


CD 
























z z 


— X 






•— « 




•— * «s 


1 — 


X O 








u 
















< <c 


H- 






> 




— > 


3 


Ld 








u 












OJ 




» • 


or cd 






<D 




•—* ui 


CD 


CD 3 




Q 
















* 






O CD 






Z 




— • o 




CD 




Ul 




X 












♦ 




o o 


3 CO 










•S- *— • 






3 


3 


3 OJ 


3 CL 
















• » 








X 




00 — 


O 


^ . — . 






«* * 


a 








OJ 






CD CD 






CO 




LO 


• 


O CD 


ft-1 




*-i * 


ft-1 CO 








W 


♦ 






• • 


« V • 










«s (D 


Ld 


• O 


It 


» 


u — 


II \ 










* 


OCJHKY-O 




tH 


CD 




-4~ X 


Z 


<3 0J 


V 




V V 


V — 




00 


ft— i 






Ul 


♦ 


CD 3 


U> Ld 




00 


z 




— X 


m- 


Ld ^ 




x 




X 




CO 


o0 




V 3 


♦ 


• • 


* • 




OJ 


ft— » 


00 


Ul 


ca 


• vO 






CO 


— 




ft— i 


0J 


»H 






LiJ CO CO 


OJ Ld CO 00 




h- 


H 


Ld X 


Ul 


O — 


00 


CO O 3 + 


in co 






-4- 


O 


CJ 


X 


Ld 


^ 'w- 


O II V 


II 


<n 


< 


h- o 


» — I 


Ul Ul ^ 


u 


O ft— < 1 i 


ft— f u 




CD 


II CD 


OJ 


+ 


+ 


» 




ru — ~ 


CO 


h- 


a: 


X 


CD 3 




D K OJ 


—ft 


O 0J u 


0J 




X 


CO X 




u 


cd 


CD 




x 


o 




Ld 


Ui 


U> 3 




— ^ 




V 


It II 


X 


It 




o 


<o 


u 


n 


ii 


u u 


tL> — U 


Ld <D 


y- 


X 


z < 


u 


u o: 


CD* 




U CD U 


CD — 


CO 


CD Ui CD 


o u 


CD 


CD 


ft — « ft— • 


O < ~ 


X CD 






Ld CJ 


ft— • 




O 


CCJ 


u ru u o co 


x cd 


X CD 












OJ 








•4- 










00 


•4- 


ID 






vO 






o 






o 








o 












ft-« 








ftH 






OJ 






OJ 








OJ 










OJ 


0J 


OJ 






OJ 



tH 

OJ 



uuu 



uuu 



47 



ITERATING BY M9VING X RIGHT 



o 

♦ o 

t w— • 
■ « 

Ui UJ 



o o 

• • 

o o 



♦ • 

Lii Ui 

* • 

' o o 

2 : Z 

< < 



o o 

• • 

o o 



in 
— c n 

(T Z 

UI u. 



% ^ 

co or 

> UJ 




_J CO 

• • 

CO CO 



ex <x 
<n go 



o o 

• « 
o o 

• • 






cu o 
















CO CO 








\ N 
















• • 






<n 


f- CO — 
















UJ UJ 






CO 


X -V h— «r-« 












X 




«*-» < — - 






cu 


uj •— 3 












■ — - 




* • 








1— *-» Q_ <D 








X 




< 




O Q 




u 


<D 


•— %. 1- h- 


o 




u 




— 


* 




z z 




— 




~ > 3 


UI 


> 


U- 


c 


X 


UJ 




<c < 




fr- 




►-« UJ CD CD 


_l 


o 




♦ 




UI 


cu 


• • 




ee 


CO 


•-« Q CO 




UJ 


1— 




CO 


1 


* 






o 


CO 




—I 


3 


JCJ -JCC 


3 X 


* 3 


-— 


_i * 


o o 


_i 


CO 




00 — O — 




’ — 


> * % C3 






X 




• • 




\ 




in • o oc «h 


X 


*-» * »-4 CO 


O 


X ^1 




h(\j y 


o o 


v~» 







CD 

-* 1“ 
— X 
UJ 
UJ »— 

o 

<D _I 
CJ _1 

z 



• »-t II + 



LJ 

zoai 

• UI > 

cr * vO 
UJ o ~ 
UI UI 



X — 



If ^ II \ 

X X X ~ 
— X 



II UJ -s II o 
X — l O X UI 



H * 



tt X 



X * O 0J *- H- X ^ C3 
CO _ 



in 



sO 

co 



— UJ _j — Ld * CO _J UuJ 

— < — X_l— X _J * • • * • 

cr, < o o i h < oin + - vex ujcocooouicnin 

^11 •COUOJI* O •OJUJXCO U OOdJUXU--OJ II X « (D 
OJ-OOJLOJ^O •OOJUJ + OJ— * • CU + + I CO — — CO »— 



x n 



II 



y ii o ii 



II UJ 



CUUQTtn-'UCDUO — CO M UJ cr> UI 



V o o 



U CO CO 



<n 



u H a) m u n Ll L ® v U UJ O 



UIU-.-- i: 0 <CU 0 U 0 < yUJUDULJOULL(.DDLGO--OfIJ-yO 



0% 


o 


rH 




in 


vO 


r-v 


00 




cu 


cu 


cu 


cu 


cu 


cu 


cu 


cu 


cu 


cu 


cu 


cu 


cu 


cu 


cu 



CJ CJ 



48 



ITERATING BY M0VING X LEFT 



NO 



— CO 

a: 2 : 

UJ Ll. 



— co 

DC Z 
Ll) Ll. 



CO oc 
•s LU 



% 

co cr 

s LU 





































Q_ 




% 


tH 












OJ 


O 




























Z> 




OJ 


O 














•r-4 


r>v 
































*-* 










X 


CO 


w 


r-4 


























X 


X 


CO 


v-4 










X 




X 


tH 




























X 




1 — rH 










UJ 


*— # 


Z> 




























0 


UJ 




z> 










X 




o_ 


CD 


























2 


K 


►r* 


a. cn 










►— « 




H- 


X 


O 




V 










u 










•— « 


•—* 




L- X 


O 




X 






> 


z> 




UJ 














u 










> 




> 


3 


UJ 




*s 




*— « 


Ui 


CD 


O 


_J 




0 










— *■ 










CD 




Ui 


CD CD 






O 




•—» 


0 




O 






LU 










1— 










21 


> — « 


0 


O 






LU 




% 


t—i 








-J 


_J 




_1 cu 


_J 


QC 














•— t 






_l 


3 




00 


■ — 


O 








■ — • 






* 




q 










X 


00 


— ~ 


O ~ 






— 




in 




• 


0 


CO 




X 




«-l 


★ 




CO 










ao 


in 




• O 


OJ 


tH 


X 






CD 


UJ 


• 


CO 


II 


* 




u 




II 


\ 














CD 


UI • 




M 


+ 






X 


2 


a oj 


X 








X 


— 




-J- 




in 


CD 




1- 


2 C 


OJ 


X 








X 


• 


UJ 


% 




x 






— • 




X 




OJ 




NO 


2 


— ^ 


X 


* UJ 






X 


vO 




UJ 


DC 


• 


vO 




— ■ 






< 




■ — 




OJ 




oj 




I s - 


LU 


cr • 


0 




■ — 


fl 


UJ 


X 


UJ 


0 






< 


0 10 


4 


vO 


<c 






sO 




*— 


n LU 


f- 


UI 0 


' — 


CO 


cn 


X 


0 




» — < 


UJ 


UJ CO 


H 


• 


CO 


U 


CO 


II 




CD 


11 


<n 


< 


X O 




► — « 1 1 1 


LU 




H 


X 


CD 




• — 


_i 


X 


OJ 




O OJ 


u 


OJ 




▼4 1 — CO 


\— 


d: 


X CD 


_1 


— 3 


f- 


OJ 




UJ 


u 








•— t 




V 


II 




II 




X 


II 




Q 




UJ 


LU LJ 


_! 








X 


h- 


2 


< 


u. 


U_ DC 


CD 




u 


CD 


u 


CD 




CO 


<r, 


LU 


CD 


»— 


X 2 


< 


L. Ll (D 






UJ 


LJ 


» — * 






0 


< 


u 0 


u 


0 


< 


X 


a 


X 


O 




• UJ 


LJ 


•— * • — • 


2: 


O CD 



U 



O -J- 
• -=r 

O OJ 
u 

Ll. <D 
Ll. O 



00 

CU 

(Vi 



CO 

cu 



10 

CO 

OJ 



SO 

CO 

OJ 



00 

CO 

OJ 



OJ 



CO 

OJ 






u u u 



49 



244 FF=FF+B(K)**2 
D9 245 K“ 1 <* L, 

245 B(K)=B(K)/SQRT(FF) 
KS«=0 



0 

• o 

w-% • 

1 *-• 
■ IT 

Ui l u 



o o 

• • 

o o 




-J _i 
• • 
LlJ Ui 



• • 

o o 
z z 
<c c 

♦ « 

o o 

• * 

o o 



— I CD 

• • 

CD CO 

• • 

or <r 

CD CD 

• • 

o o 

• • 

o o 



z 



00 
— CO 

O' z 

UJ u 
►— « — 

co cr 

•K UJ 





























%- 1- 












CD 








«— < 












































CD CD 












O 






OJ 


o 












































• • 




















tH 


rN 




























. — . 














UJ UJ 






o 






X 






m 
































X 














- — ^ 






OJ 










X 




X 






























■ — 














• • 












CD 




UI 


•— « 


Z> 
























X 






CO 














o o 




u_ 


<d 






Z 




h- 


>-4 


cl 


CD 
























— 




♦ 














z z 






1 — 






*— * 




» — 1 






*— 


Q 




X 
















CO 






UJ 














< < 




h- 








> 




— 


> 


r> 




UJ 




















♦ 


«. — 




UI 










0J 




• • 




QC 


CD 






CD 




•— « 


UJ 


CD 


CD 


_J 




o 


















<c 




1 










* 








GJ 


CD 






X 




►H 


o 




CD 






UI 














X 


* 










_J 




* 




o o 


_J 


CO 












•V 


» — « 










_J 




_J OJ 


_J 














X 














• • 




\ 








X 




00 




O 


— 


**** 










* 








*-r 


o 


V -H 


•* 






^ OJ 


X 




O O ^ 


— 








CO 




in 




•> 


O 


r^ 




X 




« — ( 


♦ 


»— 1 






H 


UI 


•V 


u 


o 






H 


♦ 






• • 


II 


X 


• 












CD 


Ui 


• 


m 


U 


» 




II 




0 






V 


_J 


o 


X 


UJ 






X 


* 


O OJ K h 


X 




e 






CD 






>— 


z 


o 


0J 


X 






X 


X 


X 










UI 




_J 








. — . 


UI 


* 


CD _J 




u 


UJ 




00 


z 






X 


• 


Ui 






X 
















X 


-J 












X 


_J 


* 


• • 




★ 


* 




OJ 


►— « 


00 




UJ 


cr 


*. 


s0 




' — 






CD 






o 


CO 


+ 




0^ 


X 






o 




— 


UJ CD CD 




UJ 


co 




*— 


n 


UJ h- 


UI 


O 


— 


•-« 


CD O 0J 


-f 


CO 


o 


• 




UJ 


X 




n 


O O LO 


CJ 


X 


UJ 




LO 


II 


X 


u 


CD 


< 


►- 


o 




» — • 


UJ 


UJ 


vO 


It 


• 


o 


u 


vO 


• 


o 


OJ 


UJ 


+ 


OJ 




• 


• 


cu 


+ 


+ 


c 


' — w* 


UJ 


. — . 




CO H- 


<r 


X 


CD 


_J 




_J 


H- OJ 


■ — 


O 0J 


U OJ 


o 


II 




It 


UJ 




X 


o 


o 




u 


cd 


CD 






X 




o 




ui 


UJ 


O _J 






» — « 




X 


II 




It 




M 


UJ 


CD 


UI 


It 


CD 




It 


It 


CD 


It 


It 


It 


u u 


CD 




u 


UJ 


CD 


h~ 


»— 


z 


c 


u 


u 


a: 


CD 


' — 


u 


CD 


u 


CD 


UI 


UI 


O 


UI 


UJ 


o o 


u 


G 

D 

F 


cd 


CD 


►— * •— t 


o 


< 


• — » 


X 


CD 


• — « 


9— 


UJ 


CJ 


• — < 


►— * 




O 


cn 


u 


o u o 


IN. 








00 




CT> 










o 








r-f 










CO 




















0J 












<1- 














in 








in 










in 














o 






V0 




OJ 








OJ 




CXJ 










OJ 








OJ 










UJ 














0J 






OJ 





uuu 



50 



-4- 



cr 



























UJ 






— CO 










































1— 






cl z 










































z 






Ul u 
















































» — * » 










































cl 






•S 










































CL 






CO CL 
















































^ UJ 






























X 








V 




UJ 






*-» •— * 






























■ — 








- — ■ 




z 
































<n 




CO 




in 




CO 




#— « 






*r-i ^ 


























vO 








00 








3 






% 


























OJ 




M 




OJ 














CU O 










































UJ 






■S, -*-* 


























CO 








CD 








X 




j— co — 


























h- 








h~ 








h- 




X 


*s ¥- 






































' — 








UJ 


~ 3 


























<D 




c 




CD 




<C 




z 




►- 


— CL 
















— v 


lu 








o 








CD 








CD 






•v \ — 
















~ 00 


Li. 






























> 3 






— * 










00 — 


■ — • 
























l~ 




* — * 


UJ CD 






vO 








— 


— X 


f- 








o 








O 








3 




•— * 


O 














X 


X CL N 


cl 
















• 




— 




a 








O 


vO 


X 




00 


X 


CL 


0 J CO 


o 








0 J4- 




CD 


_J o 




h- 




00 


— o 


nO 


CD 


CL 




t — 4 


CL 


3 


3 C ^ 


cn 








U! 








UJ 




ca 




3 




in 


m 




■% 


3 






3 < 


< u O 










♦ 


t cu 




• 




OJ 




CD 






CD UJ 0J OJ 


< 




T-l 


< 


« 


M — OJ 


— * 


hv 






CD 


II 






CD 


II 












1 - z 


H 


ii 


ii 




11 


U 




^ (J1 II 




-3- 




00 


jy vO ci 


j y oci 


CD 




’ 


X • 


* — * 


3 




Ul 


3 


■ — 


3 0A 3 




OJ 




OJ 


Li. 






vC 


u 




— 


vO 


Z 


*— < 




LU CL 






3 


3 




3 






CD 




00 




CL CO 






a 


(JA 






•— i 


II 


UJ 


f- Ul 


«— * in 




Z 


sO 






•* <r-» CTA 


N 


CD 


ii 


<n 


z r- 


LlJ CD 


z 


00 


LU 


CD 


h- 


1 


o 


►— « 


o 


O 


*—• 


*— « 


O 


* — < 


v0 


i vO O 




^ J — c o 1 — 


— 


OJ 


h- 


1— 




0J 


K- 


I — 


h- 


X 


CD 


3 — 


O) CO 




H- 


CO 


— 


•w 


CD 




II 


Q 








•— • 








»— • 




CD 


uj u; 


3 






U 


Z 




u u u u 




cd <n 


UJ CD 


U CD 


cl 


CD 


U_ CD 


CL 


CD 


3 


►— 


z 


< U CD 


CD 


*— • 


CD 


CD 


8—. 


*— • 


(D 


CD 


x cd 




o 


*— « 


o 


3 : 


CD 




O 


3 


CD 


CL 




UJ o — 


QQG-UOa 


u 


CL CL O 


CO 








in 




CO 


oa *-» 




cri in 














in 


*— 4 






V0 




O 




vO 




nO 






NO 


00 




00 


00 




0\ 










o o 






o 




CVj 




oj 




cu 




OJ OJ cu 




0J OJ 




cu 










CD CD 






CD 





u id o 



51 



P I C < 1 > J)=ALPH<7> 
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333 FORMAT ( ' 1 » ) 

338 FORMAT (20X/61A1) 



352 FORMAT (5X> * REGRESSING TS LAST STEP') 
ST0P 
END 







f 



T 
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APPENDIX II. PROGRAM OPERATING INSTRUCTIONS 



There exists two methods by which this program can be 
executed on the XDS 9300. The first and fastest method is 
utilizing a magnetic tape. The memory dump of the compila- 
tion from the digital computer has been recorded onto a 
magnetic tape. The elapsed time until execution commences 
would be approximately five minutes using this tape. The 
other method is simply reading the Fortran program from 
the computer cards. It would take roughly thirty minutes 
until the program is ready to execute using the cards. This 
time factor is very crucial if more than one set of data 
cards must be read. 

In executing the program from the magnetic tape, the 
operator may consult the publication titled "Electrical 
Engineering Computer Laboratory" by R. D. Delaura [2]. A 
brief guide in using the magnetic tape mechanism applicable 
to the object language filed on this, tape follows: 

(1) Load the magnetic tape 

(2) Unit select to position "1" 

(3) Automatic/manual switch to "MANUAL" 

(4) Depress forward drive (tape will stop at 
the LOAD POINT) 

(5) Automatic/manual switch to "AUTOMATIC" 

The magnetic tape is now ready for instructions from the 
XDS 9300. The operator must incorporate five control cards 
into the data deck. These control cards are instructions 
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to the various systems used in the execution of this program. 
The deck of cards will be as follows: 

(1) B00T 

(2) J0B 

(3) AGT 

(4) ASSIGN XI = MT1A 

(5) RERUN 

(6) (DATA DECK) 

The deck is now ready to load into the card reader and 
execution can commence. 

The following instructions will be applicable to the 
use of computer cards. The list of control cards, that 
must precede the source program printed in Appendix I, 
are given below: 

(1) B00T 

(2) $ DATA 

(3) $PATCH 

(4) 007067 001 06711 

(5) $END 

(6) J0B 

(7) AGT 

(8) FORTRAN LS , G0 

The data deck is placed at the end of the Fortran program. 

The deck of cards is now ready to load and execution can 
commence . 

In preparing the computer program using either method, 
four cards must be punched and inserted into the data cards. 
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The first three cards can be used for the purpose of verbally 
describing the data to be investigated. This information 
can be punched using any part or all of the eighty column 
card. The operator must insure in insert these three cards even 
though one or more cards may be blank. The fourth card is the 
instruction by which the computer reads the data from the 
punched cards. For illustrative purposes. The statement might 
read, " (4F7 ,2 ,3X ,3F9 . 2) . " This indicates that a total of 
seven parameters will be read from each card. The first 
parenthesis must be punched in column seven of the card. 

The actual data cards are inserted behind these four cards. 

The four cards plus the cards of data will be considered 
the data deck. 

A guide for operating the XDS 9300 computer and the A6T 
console using either a magnetic tape or cards will be given. 
First though, assume that there are M number of X points 
and N number of Y points. The value of M and N must be 
less than or equal to eighty. In addition, all observations 
in both sets constitute an L-dimensi onal vector. The number 
of L must be less than or equal to thirteen. Now, the 
instructions will be presented. 

(1) Use the "POWER ON" procedure for the XDS 9300 
and the AGT graphic systems. 

(2) Executing the program 

(a) At the XDS 9300 control panel, depress "SENSE 
2" if using cards. Do not depress if using the 
tape. Then depress "IDLE", "CLEAR" and "CLEAR 
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FLAGS" simultaneously, "RESET", "RUN", and 
"CARDS." The card reader will commence reading 
the punched cards. The line printer and 
teletypewriter (TTY) will simultaneously print 
"JOB." The computer will now start compiling 
the program. 

(b) At the time execution of the program starts, the 

TTY will print the message: "TYPE IDEV = 2* AND 

A CARRIAGE RETURN IF USING AGT2 . " The input 
light will illuminate to alert the operator. 

The operator must type the number of X's and Y's 
plus the dimension of the two subsets of data. 

In addition, the operator has the option of 
printing all WRITE statements prior to statement 
#57 in the computer program. This is done by 
typing "NPFLG = 1" if a printout is desired; 

NPFLG = 0" for no printout. As an example, 
suppose there are 15 X's and 20 Y's in 10 
dimensions. If the operator is using AGT 2 and 
requests no printout, type "M = 15, N = 20, L = 10, 
NPFLG = 0, IDEV = 2, *" and a carriage return at 
the computer TTY. If using AGT 1, simply delete 
"IDEV = 2." The card reader will begin the 
process of reading the card. 

(c) The word "NAMELIST" and a blinking cursor will 
be displayed on the bottom left corner of the 

AGT screen. This will occur when computer reaches 
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statement #57 in the program. The operator 
must type in the components of the A and B 
vectors. The values from Appendix IV may be 
used or the operator may use any two orthonormal 
unit vectors. Thirteen numbers for each vector 
must be inputted. If the sample is 10-dimensional, 
zero.es will be inserted for the last three 
components. A representation of the above example 
would look like: "A = .31628, .31628,..., .31628, 

0,0,0" and a carriage return, "B = .31628,..., 
-.31628,0,0,0" and a carriage return, and a 

carriage return. The XDS 9300 will compute the 
projections and relay the data to the AGT for 
display on the screen. The "X's" and "0's" denote 
the original "X's and Y's, respectively. 

(3) Interations and output 

(a) To change the orientation of the plane, the 
operator chooses either the "PICK X" or "PICK 

X / 

Y" function switches. The location of these 
buttons on the overlay are shown in Figure 
A-2-1 . A video message, "TEXT EDIT REQUEST," 
will appear at the bottom of the screen. The 
operator depresses the button on the lightpen 
and places the light over the desired symbol 
(point) to be moved. Upon releasing the button, 
a portionof the letter will disappear to 
signify the selected point. The user pushes 
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the two switches "COMPLETE CHOICE" and "END 
EDIT". The video message will disappear. The 
iteration will commence after one of the move- 
ment buttons (function switches) is depressed: 
"RIGHT", "LEFT", "UP", "DOWN." The 9300 
recomputes the projections from the inner product 
with the orthonormal vectors and displays the new 
results on the screen. 

(b) A printout of the new projections can be obtained 
by simultaneously depressing the "PRINT" and one 
of the movement function switches. The movement 
switch must be released first. A copy of the new 
projections on the plane will be printed at the 
line printer. If the plot of the computer 
projections is desired, then depress the "PLOT" 
function swi tch . 
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APPENDIX III . DATA LISTING 



This appendix presents a listing of the data sets used 
in the seven cases discussed in Section IV. The data has 
been arranged so that a row of numbers represents a single 
vector. For example, the data set for Case one is presented 
in three columns denoting a three-dimensional entity. The 



thirteen columns in 


Case Seven 


represents a thirteen- 


dimensional problem 


. As an aid 


for reference purposes 


following index is 


prov i ded : 






NUMBER 


NUMBER 


PAGE 


OF Y'S 


OF X'S 


CASE 1 


21 


' 0 


CASE 2 


20 


20 


CASE 3 


70 


70 


CASE 4 


40 


20 


CASE 5 


50 


25 


CASE 6 


38 


44. 


CASE 7 


60 


50 
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CASE 



1 DATA 



I YU/1) 



1 


1 


♦ 0000 


2 


1 


♦ 0000 


3 


2 


♦ 5000 


A 


2 


♦ 5000 


5 


4 


.0000 


6 


4 


.0000 


7 


1 


.0000 


8 


1 


.0000 


9 


2 


.5000 


10 


2 


.5000 


11 


4 


.0000 


12 


4 


.0000 


13 


1 


.0000 


14 


1 


• oooo 


15 


4 


.0000 


16 


4 


.0000 


17 


2 


.5000 


18 


1 


.7500 


19 


1 


.7500 


20 


3 


.2500 


21 


3 


.2500 



(1*2) Y ( I * 3 ) 



0000 


1.0000 


0000 


2.0000 


0000 


1.0000 


0000 


2.0000 


0000 


1.0000 


0000 


2.0000 


0000 


1.0000 


0000 


2.0000 


0000 


1.0000 


cooo 


2.0000 


0000 


1.0000 


0000 


2.0000 


0000 


1.0000 


0000 


2.0000 


0000 


1.0000 


0000 


2.0000 


0000 


3.3229 


5000 


2.6614 


5000 


2.6614 


5000 


2.6614 


5000 


2.6614 



Y 

1 

1 

1 

1 

1 

1 

3 

3 

3 

3 

3 

3 

2 

2 

2 

2 

2 

1 

2 

1 

2 
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CASE 2 DATA 



I 


Y < I > 1 > 


1 


3.0000 


2 


4.0000 


3 


5.0000 


4 


3.0000 


5 


4.0000 


6 


5.0000 


7 


5.0000 


8 


3.0000 


9 


3.0000 


10 


4.0000 


11 


5.0000 


12 


3.0000 


13 


3.0000 


14 


4.0000 


15 


5.0000 


16 


5.0000 


17 


5.0000 


18 


5.0000 


19 


3.0000 


20 


3.0000 



OJ 


Y ( I / 3 ) 


.0000 


5.0000 


.0000 


5.0000 


.0000 


5.0000 


.0000 


5.0000 


.0000 


5.0000 


.0000 


5.0000 


.0000 


5.0000 


.0000 


5.0000 


.0000 


7.0000 


.0000 


7.0000 


.0000 


7.0000 


.0000 


7.0000 


♦ 0000 


7.0000 


.0000 


7.0000 


.0000 


7.0000 


.0000 


7.0000 


.0000 


6.0000 


.0000 


6.0000 


.0000 


6.0000 


.0000 


6.0000 



Y 

2 

2 

2 

4 

4 

4 

3 

3 

2 

2 

2 

3 

4 

4 

4 

3 

2 

4 

4 

2 
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vl 

1 

a 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

ao 



CASE 



a 



( C 0 NT ) 



X ( J , 1 > 

7.0000 

8.0000 

9.0000 

7.0000 

8.0000 

9.0000 

9.0000 

7.0000 

7.0000 

8.0000 

9.0000 

7.0000 
7-0000 

8.0000 

9.0000 

9.0000 

9.0000 

9.0000 

7.0000 

7.0000 






X < J / 2 ) 


X ( J / 3 ) 


6.0000 


9.000 


6.0000 


9.000 


6.0000 


9.000 


8.0000 


9.000 


8.0000 


9.000 


8.0000 


9.000 


7.0000 


9.000 


7.0000 


9.000 


6*0000 


11.000 


6.0000 


11.000 


6.0000 


11.000 


7.0000 


11.000 


8.0000 


11.000 


8.0000 


11.000 


8.0000 


11.000 


7.0000 


11.000 


6.0000 


10.000 


8.0000 


10.000 


8.0000 


10.000 


6.0000 


10.000 
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CASE 3 DATA 



I 


Y ( I / 1 ) 


Y( 1/2) 


Y( 1/3) 


Y( 1/4) 


Y(I/5) 


Y( 1/6 


1 


1.180 


1.718 


1.100 


1*892 


1.617 


1.703 


2 


1.597 


1.621 


1.526 


1.358 


1.552 


1.746 


3 


1.197 


1.597 


1.540 


1.590 


1.733 


1.484 


A 


1.506 


1.561 


1.035 


1.176 


1.892 


1.899 


5 


1.588 . 


1.328 


1.639 


1 . 625 


1.754 


1.684 


6 


1.018 


1.094 


1.289 


1.646 


1.284 


1.520 


7 


1*517 


1.627 


1.317 


1*812 


1.197 


1.137 


8 


1.701 


1.452 


1.342 


1.452 


1.870 


1.917 


9 


1.124 


1.239 


1.711 


1.215 


1.546 


1.209 


10 


1.883 


1.120 


2.000 


1.072 


1.344 


1.397 


11 


1.649 


1.420 


1.998 


1.861 


1.498 


1.810 


12 


1.032 


1.365 


1.913 


1.283 


1.460 


1.825 


13 


1.929 


1.224 


1.880 


1.747 


1.405 


1.080 


14 


1.064 


1.996 


1.541 


1.837 


1.046 


1.039 


15 


1.069 


1.291 


1.979 


1.638 


1.170 


1.195 


16 


1.092 


1.226 


1 . 174 


1.806 


1.244 


1.802 


17 


1.524 


1.337 


1.919 


1*166 


1.532 


1.948 


18 


1.176 


1.690 


1 .487 


1.687 


1.422 


1.413 


19 


1.585 


1.648 


1.962 


1.225 


1.378 


1.627 


20 


1.632 


1.175 


1.615 


1.571 


1.724 


1.216 


21 


1.171 


1.659 


1.643 


1.252 


1.842 


1.958 


22 


1.684 


1.932 


1.437 


1*162 


1.378 


1.827 
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I 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 



Y < I # 1 ) 

1.208 
1.968 
1.490 
1.692 
1 * 530 
1.003 
1.344 
1.649 
1.461 
1.701 
1.745 
1.386 
1.757 
1.122 
1.289 
1.848 
1.519 
1.312 
1.777 
1 * 504 
1.489 
1.599 



CASE 3 (C8NT) 



Y( 1/2) 


Y ( I # 3 ) 


Y< 1/4) 


YU/5) 


Y( 1/6) 


1.206 


1.544 


1.267 


1.497 


1.653 


1.342 


1.531 


1.487 


1.153 


1.074 


1.133 


1.575 


1.929 


1.084 


1.715 


1.328 


1.899 


1.830 


1 . 106 


1.222 


1.909 


1.522 


1.227 


1.805 


1.239 


1.158 


1 .080 


1.096 


1.783 


1.462 


1.059 


1.792 


1.204 


1.515 


1.290 


1.484 


1.656 


1.244 


1.874 


1.487 


1.181 


1.636 


1.397 


1.924 


1.152 


1.317 


1.784 


1.409 


1.747 


1.171 


1 . 655 


1.551 


1.389 


1.368 


1.882 


1.039 


1.643 


1.236 


1.158 


1.472 


1.371 


1.398 


1 . 550 


1.861 


1.213 


1.396 


1.045 


1.647 


1.294 


1.450 


1.639 


1.797 


1.641 


1.825 


1.483 


1.173 


1.978 


1.774 


1.545 


1.693 


1.047 


1.679 


1.418 


1.915 


1.388 


1.741 


1.703 


1.528 


1.054 


1.775 


1.019 


1.415 


1.232 


1*442 


1.839 


1.072 


1.080 


1.236 


1.507 


1.236 


1.230 


1.244 


1.393 


1.287 


1.450 


1.176 


1.485 


1.616 


1.422 


1.633 
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CASE 3 ( C 8 NT ) 



I 


Y ( I / 1 ) 


Y ( 1 / 2 ) 


Y ( 1 / 3 ) 


YU / 4 ) 


Y ( 1 / 5 ) 


Y ( 1 / 6 ) 


45 


1.997 


1.119 


1.400 


1.364 


1.144 


1.194 


46 


1.239 


1.355 


1.864 


1.233 


1.052 


1.657 


47 


1.842 


1.195 


1.941 


1.176 


1.052 


1.950 


48 


1.724 


1.840 


1.602 


1.833 


1.907 


1.756 


49 


1.007 


1.528 


1.540 


1.648 


1.532 


1.372 


50 


1.952 


1.902 


1.775 


1.218 


1.930 


1.344 


51 


1.668 


1.392 


1.636 


1.575 


1.292 


1.937 


52 


1.395 


1.263 


1.838 


1.016 


1.905 


1.538 


53 


1.492 


1.989 


1.514 


1.500 


1.607 


1.803 


54 


1.629 


1.333 


1.027 


1.040 


1.741 


1.282 


55 


1.977 


1.646 


1 .145 


1.948 


1.372 


1.585 


56 


1.474 


1.988 


1.088 


1.578 


1.882 


1.557 


57 


1.904 


1.334 


1 . 104 


1.974 


1.088 


1.021 


58 


1.184 


1.770 


1.259 


1.116 


1.551 


1.654 


59 


1.029 


1.262 


1 .112 


1.501 


1.344 


1.350 


60 


1.582 


1.701 


1.598 


1.656 


1.025 


1.647 


61 


1.592 


1.496 


1.963 


1.020 


1.805 


1.703 


62 


1.296 


1.101 


1.369 


1.092 


1.789 


1.297 


63 


1.059 


1.146 


1.531 


1.058 


1.992 


1.649 


64 


1.165 


1.174 


1.485 


1.724 


1.140 


1.005 


65 


1.430 


1.907 


1.089 


1.360 


1.180 


1.028 


66 


1.746 


1.939 


1.344 


1.296 


1 .626 


1 . 260 
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CASE 3 ( C0NT ) 



I 


YU;1) 


YU/2) 


Y U / 3 ) 


YU/4) 


YU, 5) 


YU/6) 


67 


1.203 


1.601 


1.525 


1.610 


1.575 


1.151 


68 


1.431 


1.050 


1.242 


1.270 


1.146 


1.874 


69 


1.286 


1.170 


1.939 


1.011 


1.486 


1.768 


70 


1.C89 


1.497 


1.623 


1.972 


1.532 


1.876 
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J 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 



CASE 3 (C8NT ) 



X(J,1) 


X( J/2) 


X ( J/ 3 ) 


X( J/4) 


X( J,5> 


X< J/6) 


.528 


. 236 


.840 


.990 


.761 


.730 


.063 


.843 


.799 


.273 


.804 


.882 


.251 


.784 


.924 


.409 


.471 


.011 


. 466 


.538 


.499 


.784 


.751 


.413 


.293 


.315 


.990 


. 536 


.259 


.084 


.023 


.178 


.375 


.947 


.090 


.274 


.239 


.141 


.505 


.949 


.500 


.762 


i 402 


.277 


.659 


.664 


.139 


.676 


.335 


.318 


.463 


.551 


.079 


.469 


.573 


.189 


.935 


.751 


.155 


.920 


.971 


.905 


.657 


.175 


.275 


.628 


♦ 949 


.428 


.674 


.400 


.400 


.171 


.967 


.284 


.305 


.450 


.688 


.055 


.928 


.429 


.543 


.477 


.709 


. 642 


.112 


.149 


.884 


.582 


.252 


.825 


.694 


.080 


.420 


.715 


*068 


.300 


.171 


.265 


.371 


.528 


.462 


.282 


.543 


.909 


.612 


.064 


.572 


. 150 


.233 


.987 


.346 


.305 


.748 


.877 


.439 


.327 


.231 


.763 


.198 


.816 


.182 


.594 


.759 


.751 


.420 


.027 


.334 


.007 


.058 


.455 


.079 


.571 



70 



J 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 



CASE 



3 (C0NT ) 



X(J/1) 


X( J/2) 


X( J/3) 


X( J>4) 


X( J/5) 


X( J/6) 


.366 


' .942 


.747 


.010 


.371 


.237 


.067 


.288 


. 326 


.129 


.143 


.353 


.304 


.483 


.219 


.195 


.713 


.821 


1.000 


.257 


.693 


.116 


.194 


.278 


.724 


.505 


.323 


.180 


.806 


.718 


.702 


.384 


.307 


.345 


.351 


.636 


.506 


.967 


.298 


.029 


.102 


. 656 


.766 


.787 


.584 


.131 


.092 


.119 


.038 


.807 


.826 


.044 


.988 


.259 


.643 


.534 


\ 

.796 


.098 


.963 


.982 


.296 


.412 


.109 


.692 


.836 


. 105 


.337 


.672 


.070 


.882 


.448 


.305 


.398 


.659 


.515 


.768 


.378 


.627 


.476 


.091 


.688 


.552 


.897 


.989 


.477 


.892 


.514 


.633 


.048 


.520 


. 1 66 


.004 


.062 


.839 


.412 


.591 


.596 


.112 


.665 


.498 


.376 


.041 


.366 


.208 


.867 


.493 


.592 


.738 


.060 


.130 


.542 


.230 


.975 


.796 


.433 


.373 


.658 


♦ 486 


.854 


.466 


.955 


.541 


.412 


.861 


.513 


.525 


.112 


.364 


.799 


.807 


.380 


.097 



71 



J 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 



CASE 



3 (C0NT ) 



X( Jil) 


X( J/2) 


X( J/3) 


X( J#4) 


X( J,5) 


X( J/6) 


.325 


.321 


.120 


.010 


.338 


.971 


.134 


.692 


.961 


•127 


.770 


.188 


.767 


.734 


.132 


.436 


.145 


.476 


.234 


.679 


.315 


.683 


.080 


.980 


.318 


.306 


.583 


. 455 


.137 


.658 


.890 


.423 


.538 


• 126 


.138 


. 626 


.621 


.712 


.449 


.148 


.108 


.278 


.421 


.105 


.245 


.284 


.934 


.728 


.718 


.146 


.622 


.702 


.523 


.575 


.826 


♦ 635 


. 266 


.808 


.058 


.599 


.940 


♦ 189 


.933 


.178 


*540 


.429 


.636 


.978 


.629 


.380 


.292 


.874 


.780 


.202 


.703 


.072 


.150 


.824 


.981 


♦ 188 


.637 


.680 


.416 


.818 


.884 


.289 


.694 


.340 


.876 


.109 


.742 


.621 


♦ 901 


.104 


.638 


.532 


.165 


.094 


.425 


.477 


.698 


.248 


.558 


.528 


.688 


♦ 278 


.953 


.264 


.570 


.175 


.717 


.916 


.320 


.253 


.511 


.066 


.864 


.498 


.872 


.186 


.559 


.044 


.175 


.705 


.016 


.562 


.403 


.931 


.299 


.654 


.600 


.733 



72 



CASE 3 (C0NT) 



J 


X(J,1) 


X( J/2) 


X( J/3) 


X( J/4) 


X( J,5> 


X( J*6> 


67 


.463 


• 

00 

»-* 


♦ 116 


.074 


.870 


.213 


68 


.739 


.094 


.651 


.967 


.349 


.173 


69 


.50^ 


. 663 


.638 


• 712 


.513 


.822 


70 


.212 


.056 


.592 


.658 


.391 


• 566 



73 



I 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 



CASE A DATA 



Y( I# 1 ) 

• • 2500 
- 1*2000 
-.5000 
1*0000 

1.0000 
1.2000 
-.5000 
.5250 

1 .1000 
1.1000 
- 1.1000 
- 1 . 1000 
1.5000 
1.5000 
.0470 
- 1.2000 
- 1.2000 

1.0000 

1.0000 
.2000 
- . 682 C 
- . 7650 



Y ( I # 2 ) 

- 1.2000 

-.0500 

- 1.0000 

-.7000 

.5000 

.7000 

1.0000 

1.0540 

1.1000 

- 1.1000 

1.1000 

- 1.1000 

.5000 

-.5000 

- 1.5000 

.5000 

-.5000 

-.4000 

-.4000 

-.1000 

-.3590 

1.4370 



Y ( I # 3 ) 

.7000 

- 1.2000 

-.5000 

-.6000 

- 1.0000 

-.7000 

.6000 

1.1500 

.2000 

-.2000 

.2000 

-.2000 

.8000 

.7500 

-.3000 

1.0000 

1.0000 

- 1 . 4 C 00 

- 1.4000 

1.0000 

-.1320 

.6250 



Y ( I # 4 ) 

- 1.0000 

-.0500 

- 1.0000 

1.0000 

.8000 

- 1.1000 

1.2000 

.0550 

.2000 

.2000 

-.2000 

-.2000 

.6000 

-.7500 

.9400 

.5000 

-.5000 

.4000 

-.4000 

1.3000 

- 1.8380 

.8080 



74 



I 

23 

24 

25 

26 

27 

28 

29 

•■1 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 



CASE 4 ( C5NT ) 



Y( lil) 


Y ( I i 2 ) 


Y( 1*3) 


Y ( 


.5220 


-.4680 


-.2980 


-1. 


• » 7580 


-'.5810- 


1.7140 




* 1 . 9940 


.0860 


-.1110 


m f 


.7020 


-.7940 


1.6710 




.3690 


- . 8250 


1.6220 




•..8440 


-.2490 


.0820 


-1. 


-.7320 


• 0710 


-.4260 


-1. 


-.7350 


.0330 


-.0700 


-1. 


-1.9130 


♦ 3420 


.1820 




1 

• 

00 

*-*> 

00 

o 


1.4470 


-.9770 




-.3800 


1.4660 


-1.1350 


-. 


-.5510 


1.5800 


-.9870 




-1.9410 


- . 3660 


.1630 




.1130 


1.5680 


-1.0680 




.0980 


1.5600 


.3780 


-1. 


-1*9060 


.2180 


-.4960 




-1.0050 


- * 4050 


1.6320 




-1.9770 


-.2900 


.0450 


m 9 



I # 4 ) 

8480 

3730 

0550 

2390 

5890 

7860 

7740 

8580 

2230 

0940 

6260 

4580 

0220 

3450 

1680 

1410 

1690 

0510 



75 



J 

l 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 



CASE 4 ( C0NT ) 



X ( Jil ) 


X ( J , 2 > 


X ( J * 3 ) 


X ( J / 4 ) 


.1370 


.6680 


-.5700 


- . 4500 


» « 2450 


* . 7640 


- . 5690 


-.0990 


« 

• 

O 


*.8740 


.2970 


.3620 


.4360 


.5580 


.4460 


• 2870 


.1970 


.5910 


.4740 


- .5380 


.0770 


••9320 


-.2750 


• 

• 

o 


.2170 


.6450 


- » 6260 


.1710 


.0970 


.5710 


- .7950 


.1120 


.2930 


.6520 


.2490 


.2060 


- . 2080 


-.7060 


-.4710 


- . 3580 


.1990 


-.8330 


.3580 


• 1330 


.2670 


-.7540 


.3010 


.2510 


.0280 


-.9320 


-.3580 


- .0480 


- » 4860 


-.7850 


.1990 


.2100 


. 2540 


.6090 


.2680 


-.5980 


.1630 


-.8880 


.2430 


.2220 


.4100 


.3490 


.3970 


.3460 


.2020 


-.7230 


- .5960 


.1500 


* . 3550 


.5350 


-.6690 


.3320 


.2280 


.5960 


.4750 


-.5370 



76 



1 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 



CASE 5 DATA 



( 1 / 1 ) 


Y ( 


1 * 2 ) 


Y < 1 * 3 ) 


Y ( I . 4 ) 


*1000 


3 . 


5000 


1.4000 


♦ 2000 


♦ 9000 


3 . 


0000 


1.4000 


.2000 


• 7000 


3 . 


2000 


1.3000 


.2000 


• 6000 


3 . 


1000 


1.5000 


.2000 


• 0000 


3 . 


6000 


1.4000 


.2000 


*4000 


3 . 


9000 


1.7000 


.4000 


.6000 


3 . 


4000 


1.4000 


.3000 


.0000 


3 . 


4000 


1.5000 


.2000 


.4000 


2 . 


9000 


1.4000 


.2000 


.9000 


3 . 


1000 


1.5000 


. 1000 


.4000 


3 . 


7000 


1.5000 


.2000 


.8000 


3 . 


400 C 


1.6000 


.2000 


.8000 


3 . 


0000 


1.4000 


• 100 Q 


.3000 


3 . 


0000 


1.1000 


.1000 


.8000 


4 . 


0000 


1.2000 


.2000 


.7000 


4 . 


4000 


1.5000 


.4000 


.4000 


3 . 


9000 


1.3000 


.4000 


. 1000 


3 . 


5000 


1.4000 


.3000 


.7000 


3 . 


8000 


1.7000 


♦ 3000 


. 1000 


3 . 


8000 


1.5000 


.3000 


.4000 


3 . 


4000 


1.7000 


.2000 


.1000 


3 . 


7000 


1.5000 


.4000 



77 



I 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 



CASE 5 < C 0 NT ) 



Y ( I j 1 ) 


Y < 1 / 2 ) 


4 « 6000 


3.6000 


5.1000 


3.3000 


4.8000 


3.4000 


5.9000 


3.2000 


6.1000 


2.8000 


6.3000 


2.5000 


6.1000 


2.8000 


6.4000 


2.9000 


7.0000 


3.2000 


6.4000 


3.2000 


6.9000 


3.1000 


5.5000 


2.3000 


6.5000 


2.8000 


5.7000 


2.8000 


6.3000 


3.3000 


4.9000 


2.4000 


6.6000 


2.9000 


5.2000 


2.7000 


5.0000 


2.0000 


5.9000 


3.0000 


6.0000 


2.2000 


6.1000 


2.9000 



Y { I # 3 ) 


Y ( I # 4 ) 


1.0000 


.2000 


1.7000 


.5000 


1.9000 


.2000 


4.8000 


1.8000 


4.0000 


1.3000 


4.9000 


1.5000 


4.7000 


1.2000 


4.3000 


1.3000 


4.7000 


1.4000 


4.5000 


1.5000 


4. 9000 


1.5000 


4.0000 


1.3000 


4.6000 


1.5000 


4.5000 


1.3000 


4.7000 


1.6000 


3.3000 


1.0000 


4.6000 


1.3000 


3,9000 


1.4000 


3.5000 


1.0000 


4.2000 


1.5000 


4.0000 


1.0000 


4.7000 


1.4000 



78 



CASE 5 



I 


Y(lil) 


45 


5.6000 


46 


6.7000 


47 


5.6000 


48 


5.8000 


49 


6.2000 


50 


5.6000 



(C0NT) 

Y ( I / 2 ) Y 

2.9000 3 

3.1000 4 

3.0000 4 

2.7000 4 

2.2000 4 

2.5000 3 



1/3) 


Y( 1/4) 


6000 


1.3000 


4000 


1.4000 


5000 


1.5000 


1000 


1.0000 


5000 


1.5000 


9000 


1.1000 



79 



J 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 



CASE 5 ( C 0 NT ) 



< J / 1 > 


X ( Jt 2 ) 


X ( J / 3 ) 


x < 


.5000 


3.2000 


5.1000 


2 . 


.4000 


2.7000 


5.3000 


1 . 


.8000 


3.0000 


5.5000 


2 . 


.7000 


2.5000 


5.0000 


2 . 


.8000 


2.8000 


5.1000 


2 . 


.4000 


3.2000 


5.3000 


2 . 


.5000 


3.0000 


5.5000 


1 . 


.7000 


3.8000 


6.7000 


2 . 


.7000 


2.6000 


6.9000 


2 . 


.0000 


2.2000 


5.0000 


1 . 


.3000 


3.3000 


6.0000 


2 . 


.8000 


2*7000 


5.1000 


1 . 


.1000 


3.0000 


5.9000 


2 , 


♦ 3000 


2.9000 


5.6000 


1 . 


.5000 


3.0000 


5.8000 


2 . 


.6000 


3.0000 


6.6000 


2 . 


.9000 


2.5000 


4.5000 


1 . 


.3000 


2.9000 


6.3000 


1 . 


’• 7000 


2.5000 


5.8000 


1 . 


.2000 


3. 6000 


6.1000 


2 . 


1.9000 


3.2000 


5.7000 


2 . 


.6000 


2 .8000 


4,9000 


2 . 



J#4) 

0000 

9000 

1000 

0000 

4000 

3000 

8000 

2000 

3000 

5000 

5000 

9000 

1000 

8000 

2000 

1000 

7000 

8000 

8000 

5000 

3000 

0000 



80 





CASE 5 


(C0NT) 






J 


X< J, 1 ) 


X( J»2) 


X( J,3) 


X(Ji4) 



23 


7.7000 


2.8000 


6.7000 


2.0000 


24 


6.3000 


2.7000 


4.9000 


1.8000 


25 


6.7000 


3.3000 


5. 7000 


2.1000 



81 



I 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 



CASE 6 OATA 



Yd/1) Y ( I / 2 ) Y< I > 3 } Y(I 



160 . 


5.1 


10.0 


28 


155 . 


5.2 


8.0 


27 


130 . 


5.0 


10.0 


26 


161 . 


5.0 


10.0 


27 


135 . 


5.0 


12.0 


27 


165 . 


5*0 


11.0 


23 


150 . 


5.0 


9.0 


29 


148 . 


4.8 


8.0 


26 


150 . 


4.5 


7.0 


31 


120 . 


4.0 


6*0 


33 


120 . 


5.1 


8.0 


32 


190 . 


4.2 


8.0 


30 


100 . 


4*4 


9.0 


35 


150 . 


4.0 


7.0 


29 


90 . 


4.6 


8.0 


30 


120 . 


4.7 


7.0 


35 


200 . 


4.3 


9.0 


30 


120 . 


4.1 


8.0 


28 


105 . 


5.0 


7.0 


27 


210 . 


5.2 


8.0 


26 


90 . 


4*0 


10.0 


25 


no . 


5.2 


11.0 


25 



4 ) 


Y ( I / 5 ) 


Y ( I / 6 ) 


Y ( I / 7 ) 


0 


70 , 


450 . 


4.10 


0 


85 , 


400 . 


3.12 


0 


75 . 


560 . 


2.96 


0 . 


70 . 


665 . 


4 * 04 


0 


88 . 


570 . 


3.05 


0 


95 . 


675 . 


3.72 


0 


90 . 


580 . 


3.83 


0 


85 . 


390 . 


3.75 


0 


60 . 


435 . 


3.05 


0 


55 . 


440 . 


2.78 


0 


56 . 


650 » 


2.89 


0 


55 . 


640 . 


3.46 


0 


48 . 


430 . 


2.98 


0 


65 . 


650 . 


2.84 


0 


70 . 


655 . 


2.96 


0 


67 . 


645 . 


2.78 


0 


62 . 


660 . 


3.20 


0 


63 . 


530 . 


3.05 


0 


64 . 


435 . 


3.55 


0 


67 . 


440 . 


3.10 


0 


68 . 


430 . 


3.41 


0 


60 . 


530 . 


2.97 



Y ( 1 / 8 ) 

21.0 

17.0 

19.0 

21.0 

23.0 

24.0 

25.0 

28.0 
20*0 

18.0 

23.0 

27.0 

26.0 

20.0 

17.0 

20.0 

23.0 

19.0 

27.0 

20.0 

21.0 

25.0 



82 



I 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 



CASE 6 ( C 0 NT > 



Y< J,l> 


Y< I 


. 2 ) 


Y( 1 / 3 ) 


Y(I 


/ 4 ) 


Y ( 1 / 5 ) 


Y < 1 / 6 ) 


Yd 


/ 7 ) 


Yd / 8 ) 


100 . 


4 


.3 


9.0 


25 


.0 


70 . 


440 . 


3 . 


62 


17.0 


70 . 


4 


.5 


8.0 


23 


.0 


64 . 


450 . 


3 . 


40 


19.0 


100 . 


4 


.8 


9.0 


27 


.0 


65 . 


355 . 


3 . 


10 


24.0 


130 . 


5 


.2 


9.0 


25 


.0 


70 . 


380 . 


2 . 


80 


29.0 


90 . 


4 


.5 


11.0 


37 


.0 


74 , 


350 . 


2 . 


95 


31.0 


80 . 


. 4 


.6 


10.0 


32 


.0 


78 . 


450 . 


2 . 


12 


32.0 


95 . 


4 


.9 


10.0 


25 


.0 


82 . 


260 . 


3 . 


11 


23.0 


70 . 


4 


.4 


12.0 


30 


.0 


85 . 


262 . 


3 . 


22 


24.0 


80 . 


4 


.2 


13.0 


36 


.0 


69 . 


260 . 


3 . 


35 


26.0 


95 . 


5 


.1 


15.0 


31 


.0 


70 . 


265 . 


3 . 


40 


29.0 


100 . 


4 


.6 


11.0 


24 


.0 


76 . 


270 . 


3 . 


54 


u> 

o 

• 

o 


95 . 


4 


.8 


10.0 


27 


.0 


74 . 


355 . 


3 . 


63 


31.0 


85 . 


4 


.7 


12.0 


25 


.0 


V 

73 , 


360 . 


3 . 


17 


ro 

VO 

• 

o 


70 . 


4 


.8 


11.0 


26 


.0 


78 . 


365 . 


3 . 


18 


30.0 


80 . 


5 


.4 


10.0 


21 


.0 


80 . 


370 . 


2 . 


90 


32.0 


85 . 


5 


.5 


13.0 


33 


.0 


81 * 


355 . 


3 . 


12 


35.0 



83 



J 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 



CASE 6 ( C 0 NT ) 



X(J,1> 


X ( J * 2 > 


X ( J / 3 ) 


X ( Ji 4 ) 


X ( J / 5 ) 


X ( J * 6 ) 


X ( J / 7 ) 


X ( Ji 8 ) 


200 # 


3.4 


10.0 


24.0 


98 . 


210 . 


4.10 


40.0 


260 # 


3.1 


8.0 


21.0 


110 . 


220 . 


3.89 


36.0 


195 # 


3.0 


9.0 


20.0 


105 . 


130 . 


2.98 


44.0 


195 # 


3.2 


9.0 


19.0 


110 . 


010 . 


3.10 


36.0 


220 # 


3.3 


10.0 


24.0 


95 . 


205 . 


3.54 . 


29.0 


220 . 


3.0 


8.0 


25.0 


90 . 


210 . 


3.78 


34.0 


190 . 


3.4 


9.0 


26.0 


96 . 


070 . 


3.36 


30.0 


285 . 


3.0 


11.0 


19.0 


100 . 


990 . 


4.23 


31.0 


300 . 


3.0 


9.0 


o 

• 

o 

CVJ 


1 02 . 


120 . 


3.75 


40.0 


225 . 


3.0 


10.0 


22.0 


1 05 » 


985 . 


3.80 


25.0 


260 . 


3.4 


8.0 


22.0 


97 . 


090 . 


2.94 


38.0 


280 . 


3.0 


8.0 


o 

• 

o 

OJ 


112 . 


200 . 


2.70 


41.0 


300 . 


3.4 


10.0 


20.0 


108 . 


835 . 


3.20 


45.0 


310 . 


3.0 


11.0 


19.0 


106 . 


055 . 


3.45 


37.0 


290 . 


3.1 


12.0 


26.0 


94 . 


240 . 


3.50 


41.0 


260 . 


3.0 


8.0 


22 . 0 


98 , 


015 . 


4.60 


29.0 


290 . 


3.3 


9.0 


25.0 


100 . 


010 . 


3*56 


32.0 


160 . 


3*1 


11.0 


20.0 


79 . 


170 . 


4.26 


34.0 


240 . 


3.5 


11.0 


20.0 


88 . 


990 . 


4.40 


42.0 


195 . 


3.1 


8.0 


21.0 


81 . 


975 . 


3.58 


44.0 


290 . 


3.4 


10.0 


19.0 


94 . 


860 . 


3.40 


40.0 


210 . 


3.5 


9.0 


22.0 


96 . 


950 . 


2.70 


47.0 



84 



J 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 



CASE 6 ( C 0 NT ) 



X ( J , l > 


X ( J > 2 ) 


X( J / 3 ) 


X ( J / 4 ) 


X ( J / 5 ) 


X ( J / 6 ) 


X ( J > 7 ) 


X ( J / 8 ) 


180 . 


3.0 


11.0 


22.0 


97 . 


990 . 


3.83 


28.0 


205 . 


2.9 


11.0 


23 * 0 


90 . 


805 . 


3.25 


38.0 


215 . 


3 . 4 . 


8.0 


o 

• 

tH 

OJ 


100 . 


700 . 


4.05 


32.0 


270 . 


3.1 


8.0 


o 

• 

o 

OJ 


111 . 


170 . 


4.60 


31.0 


290 . 


3.0 


9.0 


23.0 


102 . 


350 . 


4.42 


33.0 


320 . 


3.2 


10.0 


19.0 


87 . 


160 . 


3.64 


38.0 


210 . 


3.0 


9.0 


18.0 


112 . 


010 . 


3.81 


45.0 


210 . 


3*0 


8.0 


21.0 


95 . 


190 . 


4.00 


40.0 


185 . 


3.4 


9.0 


25.0 


96 . 


055 . 


4.18 


41.0 


200 . 


3.2 


8.0 


26.0 


98 . 


980 . 


3.46 


38.0 


170 . 


2.9 


9.0 


20.0 


95 . 


095 . 


4.15 


40.0 


140 . 


3.0 


9.0 


20*0 


98 . 


990 . 


4.25 


32.0 


no. 


4.9 


9.0 


22.0 


130 . 


220 . 


4.51 


19.0 


100 . 


5.6 


8.0 


19.0 


128 . 


216 . 


2.60 


14.0 


95 . 


4.9 


8.0 


24.0 


124 . 


218 . 


2.45 


10.0 


65 . 


6*2 


9.0 


30.0 


134 . 


200 . 


2.61 


12.0 


55 . 


5.0 


10.0 


27.0 


128 . 


205 . 


2.62 


11.0 


70 . 


5*3 


7.0 


28.0 


118 . 


204 . 


2.35 


12.0 


85 . 


4.9 


11.0 


19.0 


117 . 


206 . 


2.90 


14.0 


115 . 


5.0 


10.0 


21.0 


122 . 


198 . 


3.25 


17.0 


95 . 


4.8 


8.0 


27.0 


114 . 


228 . 


2.65 


11.0 


120 . 


6*1 


9.0 


24.0 


120 . 


244 . 


2.82 


18.0 
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I 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 



Y 1 3 

39 

38 

42 

42 

41 

43 

50 

44 

28 

47 

39 

37 

39 

50 

48 

35 

48 

38 

39 

49 

38 

38 



CASE 7 DATA 



Y 1 


Y2 


Y3 


Y4 


Y5 


Y6 


Y7 


Y8 


Y9 


Y 1 0 


Y 1 1 


Y 12 


1 


365 


152 


90 


4 


2 


4 


3 


1 


1 


5 


250 


1 


517 


126 


70 


1 


2 


3 


2 


2 


2 


2 


326 


1 


343 


180 


86 


1 


1 


3 


2 


2 


3 


2 


278 


1 


432 


110 


65 


4 


1 


3 


2 


4 


4 


2 


278 


1 


376 


115 


80 


1 


1 


3 


2 


1 


3 


2 


278 


1 


410 


110 


70 


1 


1 


3 


2 


2 


3 


2 


255 


1 


408 


105 


65 


1 


1 


3 


2 


2 


3 


2 


278 


1 


410 


150 


100 


1 


2 


4 


2 


2 


4 


2 


278 


1 


315 


118 


76 


1 


2 


3 


2 


2 


4 


5 


315 


1 


546 


150 


70 


4 


1 


3 


2 


2 


3 


1 


212 


1 


437 


120 


70 


4 


2 


3 


2 


2 


4 


2 


372 


1 


506 


106 


72 


4 


2 


3 


2 


2 


3 


1 


1 45 


1 


338 


110 


70 


2 


2 


3 


2 


1 


3 


1 


278 


1 


446 


120 


70 


1 


1 


2 


2 


1 


4 


1 


278 


1 


426 


135 


80 


1 


1 


3 


2 


1 


4 


1 


278 


1 


410 


140 


85 


1 


1 


3 


2 


2 


4 


1 


250 


1 


483 


120 


70 


1 


2 


3 


2 


2 


4 


1 


278 


3 


377 


114 


70 


2 


2 


2 


2 


2 


4 


5 


350 


3 


410 


130 


85 


2 


1 


3 


2 


2 


4 


2 


264 


3 


414 


120 


80 


4 


2 


1 


2 


1 


3 


5 


246 


1 


369 


140 


90 


4 


1 


3 


2 


1 


4 


1 


278 


1 


335 


132 


82 


3 


1 


3 


1 


1 


4 


1 


278 



86 



Y1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 



Y 1 3 

40 

49 

53 

45 

38 

47 

36 

41 

43 

38 

40 

47 

43 

33 

43 

51 

40 

45 

46 

41 

45 

52 



Y2 


Y3 


Y4 


CASE 

Y5 


7 

Y6 


(C0NT) 
Y7 Y8 


Y9 


Y10 


Y 1 1 


Y 12 


41 C 


140 


80 


3 


1 


3 


1 


2 


4 


2 


268 


422 


152 


92 


1 


2 


2 


2 


2 


3 


2 


278 


517 


no 


80 


1 


2 


3 


2 


2 


4 


5 


250 


250 


130 


70 


1 


2 


3 


2 


2 


4 


2 


315 


466 


114 


74 


1 


1 


2 


2 


2 


4 


3 


250 


366 


120 


90 


1 


1 


2 


2 


2 


4 


3 


.278 


271 


110 


70 


1 


1 


3 


2 


2 


4 


5 


295 


467 


130 


70 


4 


1 


3 


2 


2 


4 


5 


278 


372 


128 


92 


4 


1 


3 


2 


2 


4 


5 


353 


300 


140 


75 


4 


2 


3 


2 


2 


4 


5 


200 


348 


110 


80 


4 


1 


2 


1 


2 


4 


2 


335 


288 


130 


80 


1 


1 


4 


2 


2 


4 


2 


264 


477 


115 


80 


1 


1 


2 


2 


2 


4 


2 


245 


406 


140 


no 


2 


2 


2 


2 


2 


4 . 


5 


300 


489 


no 


64 


1 


2 


2 


2 


2 


3 


5 


278 


406 


130 


80 


4 


2 


4 


2 


2 


3 


5 


278 


367 


135 


85 


1 


2 


3 


2 


1 


3 


3 


350 


444 


125 


80 


1 


2 


3 


2 


1 


3 


1 


302 


402 


122 


70 


1 


1 


1 


2 


2 


3 


1 


278 


462 


120 


78 


1 


2 


3 


1 


2 


4 


1 


278 


511 


140 


95 


4 


1 


2 


2 


1 


3 


1 


350 


410 


140 


90 


4 


1 


2 


2 


1 


3 


2 


245 



87 



Yl 

1 

1 

1 

1 

t 

2 

2 

1 

1 

1 

1 

1 

1 

1 

1 

3 



Y13 

40 

37 

40 

42 

23 

42 

47 

47 

46 

41 

43 

43 

44 

43 

39 

52 



CASE 7 (CONT) 



Y2 


Y3 


Y4 


Y5 


Y6 


Y7 


Y8 


Y9 


Y 1 0 


YU 


Y 12 


446 


90 


60 


1 


1 


4 


1 


1 


3 


2 


202 


468 


110 


70 


1 


2 


3 


2 


1 


4 


1 


314 


410 


130 


80 


1 


2 


3 


2 


2 


4 


5 


351 


410 


140 


88 


1 


2 


3 


2 


2 


4 


1 


242 


355 


120 


70 


1 


1 


3 


2 


2 


4 


3 


366 


394 


116 


80 


1 


2 


1 


2 


2 


4 


2 


283 


312 


100 


60 


1 


2 


1 


2 


2 


4 


2 


310 


383 


100 


60 


1 


2 


1 


2 


2 


4 


2 


310 


489 


140 


80 


1 


2 


2 


2 


1 


4 


2 


278 


410 


130 


82 


1 


2 


2 


2 


l 


4 


2 


250 


404 


140 


98 


1 


2 


2 


2 


l 


4 


2 


278 


378 


140 


90 


1 


2 


3 


l 


2 


3 


1 


267 


369 


130 


70 


2 


1 


4 


2 


1 


3 


1 


278 


388 


133 


90 


■ 'i 

3 


1 


3 


2 


3 


4 


5 


288 


483 


110 


80 


4 


1 


2 


2 


2 


4 


2 


277 


461 


140 


90 


i 


2 


4 


1 


1 


2 


1 


330 



88 



XI 

1 

1 

1 

2 

l 

3 

l 

l 

3 

2 

1 

1 

1 

2 

1 

1 

1 

2 

1 

2 

2 

1 



X 1 3 

67 

43 

37 

40 

37 

52 

37 

40 

40 

42 

32 

40 

34 

41 

40 

23 

42 

41 

36 

40 

44 

43 



CASE 7 (C0NT) 



X2 


X3 


X4 


X5 


X6 


X7 


X8 


X9 


X10 


xu 


X 12 


404 


160 


78 


4 


2 


1 


2 


1 


4 


l 


200 


418 


118 


76 


1 


1 


3 


3 


1 


1 


2 


298 


288 


128 


86 


1 


2 


4 


3 


1 


1 


2 


209 


349 


170 


no 


3 


2 


1 


1 


4 


1 


5 


343 


396. 


118 


72 


3 


1 


2 


3 


1 


1 


2 


221 


411 


180 


120 


1 


2 


4 


1 


l 


1 


2 


196 


443 


160 


90 


4 


2 


3 


3 


1 


1 


2 


271 


479 


140 


90 


1 


2 


2 


1 


1 


1 


1 


260 


388 


140 


90 


1 


2 


3 


1 


1 


1 


5 


281 


445 


130 


88 


4 


2 


4 


3 


1 


1 


1 


215 


389 


110 


78 


3 


2 


4 


3 


1 


1 


2 


249 


406 


114 


70 


2 


2 


3 


1 


1 


1 


5 


248 


431 


100 


60 


4 


2 


3 


1 


1 


1 


1 


224 


413 


144 


86 


4 


2 


4 


3. 


1 


1 


1 


247 


411 


110 


50 


4 


1 


4 


1 


1 


1 


2 


240 


406 


no 


76 


1 


1 


3 


3 


1 


4 


2 


248 


444 


120 


80 


1 


2 


1 


3 


1 


1 


1 


151 


375 


120 


70 


1 


2 


3 


3 


1 


1 


3 


248 


383 


120 


80 


2 


1 


3 


3 


1 


1 


1 


256 


439 


114 


80 


1 


2 


4 


3 


1 


1 


2 


248 


489 


104 


70 


1 


2 


3 


3 


1 


3 


1 


221 


387 


118 


90 


1 


1 


4 


3 


1 


1 


2 


279 



89 



XI 

3 

1 

3 

2 

3 

1 

2 

1 

1 

1 

1 

1 

2 

1 

3 

1 

1 

1 

1 

1 

1 

1 



X2 


X3 


X4 


CASE 

X5 


7 

X6 


( C0NT ) 
X7 X8 


X9 


X10 


Xll 


X 12 


X 1 3 


351 


120 


90 


4 


2 


3 


3 


2 


5 


5 


258 


36 


350 


110 


66 


1 


1 


3 


3 


1 


1 


3 


248 


41 


536 


124 


84 


4 


2 


2 


3 


1 


1 


5 


282 


33 


432 


110 


70 


1 


1 


4 


1 


1 


2 


1 


240 


55 


393 


150 


90 


3 


1 


3 


1 


1 


1 


5 


245 


46 


369 


132 


80 


1 


2 


2 


3 


1 


1 


1 


237 


40 


418 


120 


80 


3 


2 


3 


3 


1 


3 


2 


236 


40 


383 


126 


84 


1 


1 


4 


3 


1 


0 


1 


221 


36 


325 


135 


80 


1 


2 


4 


3 


1 


5 


5 


175 


30 


395 


120 


78 


4 


1 


4 


1 


1 


1 


1 


178 


40 


375 


117 


75 


1 


1 


4 


3 


1 


1 


1 


229 


39 


432 


142 


80 


4 


2 


4 


3 


1 


1 


1 


235 


43 


267 


120 


84 


4 


2 


4 


3 


1 


1 


1 


213 


40 


276 


158 


100 


1 


2 


4 


3 


1 


1 


2 


167 


41 


489 


188 


70 


1 


2 


3 


3 


. 1 


1 


3 


248 


44 


369 


120 


84 


1 


1 


4 


3 


1 


1 


1 


294 


43 


398 


150 


90 


4 


2 


4 


3 


1 


1 


1 


223 


42 

• 


442 


118 


70 


3 


2 


2 


3 


1 


1 


5 


240 


44 


394 


130 


80 


1 


2 


3 


3 


1 


1 


5 


288 


31 


444 


120 


78 


3 


2 


2 


3 


1 


2 


1 


256 


42 


389 


140 


95 


1 


1 


4 


3 


2 


2 


5 


273 


44 


359 


130 


100 


1 


1 


2 


3 


1 


1 


5 


195 


35 



90 



J 


XI 


X2 


X3 


X4 


CASE 

X5 


7 

X6 


(C8NT) 
X7 X8 


X9 


X10 


Xll 


X 12 


XI 3 


45 


2 


232 


142 


100 


1 


1 


4 


1 


1 


1 


5 


325 


55 


46 


1 


379 


106 


72 


4 


2 


4 


3 


1 


1 


1 


264 


34 


47 


1 


378 


140 


100 


4 


1 


4 


1 


1 


1 


1 


309 


39 


48 


1 


381 


150 


90 


4 


1 


4 


1 


1 


1 


1 


282 


54 


49 


2 


364 


128 


82 


4 


1 


4 


3 


X 


1 


5 


250 


33 


50 


2 


437 


140 


90 


4 


2 


4 


3 


1 


1 


5 


248 


36 
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APPENDIX IV. TABLE OF A AND B VECTORS 



The following table is the values of the various components 
of the A and B vectors used in the seven examples. Only six 
of the thirteeen components have been listed. To use the 
table, repeat the pattern of the values until the number of 
components equals the dimension of the problem. For the 
remaining components, use a zero value until one has filled 
the thirteen components. 

TABLE OF A AND B VALUES 
DIMENSION COMPONENTS 







1 


2 


3 


4 


5 


6 


3 


A 


. 50000 


.7071 1 


.50000 


.00000 


.00000 


.00000 




B 


.50000 


- .7071 1 


.50000 


.00000 


.00000 


.00000 


4 


A 


. 50000 


. 50000 


.50000 


.50000 


.00000 


.00000 




B 


. 50000 


- .50000 


. 50000 


-.50000 


.00000 


.00000 


5 


A 


.40825 


. 50000 


.40825 


.50000 


.40825 


.00000 




B 


.40825 


- . 50000 


.40825 


-.50000 


.40825 


.00000 


6 


A 


.40825 


.40825 


.40825 


.40825 


.40825 


.40825 




B 


.40825 


-.40825 


.40825 


-.40825 


.40825 


- .40825 


7 


A 


.35355 


.40825 


.35355 


.40825 


.35355 


.40825 




B 


.35355 


-.40825 


.35355 


-.40825 


. 35355 


- .40825 


8 


A 


.35355 


.35355 


.35355 


.35355 


.35355 


.35355 


B 


.35355 


-.35355 


.35355 


-.35355 


.35355 


- .35355 


9 


A 


.31628 


.35355 


.31628 


.35355 


.31628 


.35355 




B 


.31628 


-.35355 


.31628 


-.35355 


.31628 


- .35355 


10 


A 


.31628 


. 31 628 


. 31628 


.31628 


. 31628 


.31628 


B 


.31628 


-.31628 


.31628 


-.31628 


.31628 


- .31628 


1 1 


A 


.28867 


.31628 


.28867 


.31628 


.28867 


.31628 


B 


.28867 


-.31628 


.28867 


- .31628 


.28867 


-. 31 628 


12 


A 


.28867 


.28867 


.28867 


.28867 


.28867 


.28867 


B 


.28867 


-.28867 


.28867 


-.28867 


.28867 


- .28867 


13 


A 


. 26726 


. 28867 


.26726 


.28867 


.26726 


.28867 


B 


.26726 


-.28867 


.26726 


- .28867 


.26726 


-.28867 
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APPENDIX V. TWO DIMENSIONAL PLOTS 



The following figures represent the plots of the computer 
printout. The values of the A and B vectors are printed below 
the plot if a discriminant function was determined. The 
components for the initial projections can be extracted from 
Appendix IV. A polaroid print is also included for each case. 
This allows for a comparison between the AGT display and the 
.computer plot to be obtained. 
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Figure IV-1 -1 , CASE 1 , STEP 1 
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Figure IV-1-2. CASE 1 , STEP 2 
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Figure IV-1-4. CASE 1 , STEP 4 
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Figure IV-1-5. CASE 1 , STEP 5 




Figure IV-1-5P. CASE 1, STEP 5 
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CD * 




Figure I V-2-1 . CASE 2 , STEP 1 



100 



* Q> 




Figure IV-2-2. 


CASE 2, STEP 2 


A ( 1 ) =-.0418 


B ( 1 ) = .7059 


A ( 2 ) = .9983 


B ( 2 ) = .0591 


A( 3) =-.0418 


B ( 3 ) = .7059 
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Figure IV-2-3. 
A ( 1 ) =-.3240 
A(2) =-.0005 
A ( 3 ) =-.9461 



CASE 2, STEP 3 
B(l) =-.0248 
B ( 2) = .9997 
B ( 3) = .0079 
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« 0 0 




Figure IV-2-4. 
A ( 1 ) = .0081 
A(2) =-.9998 
A( 3) =-.0161 



CASE 2, STEP 5 
B ( 1 ) =-.8980 
B ( 2 ) =-.0002 
B ( 3 ) =-.4400 
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Figure I V - 2 - 5 . 


CASE 2, STEP 6 


A ( 1 ) =-.9122 


B ( 1 ) = .0887 


A ( 2 ) = .2671 


B ( 2 ) = .8690 


A ( 3 ) =-.3107 


B ( 3) = .4867 
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Figure IV-2-5P, CASE 2, 



STEP 6 



1 05 




Figure IV-3-1 . CASE 3, STEP 1 
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Figure IV-3-2. CASE 3, STEP 2 
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Figure I V - 3 - 3 . 



CASE 3, STEP 3 




Figure IV-3-4. CASE 3, STEP 4 
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Figure I V - 3 - 5 



CASE 3, STEP 5 




Figure IV-3-6P. CASE 3, STEP 7 
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a> * 




Figure IV-4-1. CASE 4, RUN 1, STEP 1 
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Figure I V - 4 - 2 . CASE 4, RUN 1, STEP 2 
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Figure IV-4-3. CASE 4, RUN 1, STEP 4 
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Figure IV-4-4. CASE 4, RUN 1, STEP 5 
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Figure IV-4-4P. CASE 4, RUN 1, STEP 
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Figure I V -4-6 . CASE 4, RUN 2, STEP 1 
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Figure IV-4-7. CASE 4, RUN 2, STEP 2 
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Figure I V - 4 
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Figure IV-4-9. CASE 4, RUN 2, STEP 4 
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Figure I V - 4 - 1 0 . 
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Figure IV-4-10P. CASE 4, RUN 2, STEP 6 
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Figure I V - 5 - 1 . CASE 5, STEP 1 
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Figure IV-5-2 
A ( 1 ) =-.7748 
A ( 2) = .1863 
A ( 3 ) =-.1346 
A(4) = .5890 



CASE 5, STEP 2 
B ( 1 ) =-.2831 
B ( 2) =-.6366 
B ( 3) = .7173 
B ( 4 ) =-.0072 
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Fi gure 


I V - 5 - 3 . 
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Fi gure 
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Figure IV-5-4P, CASE 5, STEP 5 
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Figure IV-5-5. 
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Figure 


I V - 5 - 6 . 
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Figure IV-6-1. CASE 6, RUN 1, STEP 1 
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Figure IV-6-2. CASE 6, RUN 1, STEP 2 
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Figure IV-6-3. CASE 6, RUN 2, STEP 4 
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Figure IV-6-4. CASE 6, RUN 2, STEP 1 
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Figure IV-6-4P. 



CASE 6, RUN 2, STEP 1 
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Figure IV-6-5. CASE 6, RUN 2, STEP 2 
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Figure IV-6-6. CASE 
A ( 1 ) =-.9019 
A ( 2 ) = .0130 
A ( 3 ) = .0246 
A ( 4 ) = .0703 
A ( 5 ) =-.1947 
A ( 6 ) =-.3781 



6, RUN 2, STEP 12 
B ( 1 ) =-.3386 
B ( 2 ) = .0193 
B ( 3) = .0578 
B ( 4 ) = .0842 
B ( 5) =-.1762 
B( 6) = .9184 
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Figure IV -6-7. 


CASE 6, 


RUN 2, 


A ( 1 ) =-.6549 


B ( 1 ) 


=-.351 1 


A ( 2 ) = .0030 


B( 2) 


= .0016 


A ( 3) =-.0039 


B ( 3) 


= .0134 


A ( 4 ) = .0252 


B ( 4 ) 


= .0384 


A ( 5 ) =-.0953 


B ( 5 ) 


=-.1795 


A( 6 ) =-.7492 


B ( 6 ) 


= .9181 



139 




Figure IV-7-1. CASE 7, RUN 1, STEP 1 
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Figure IV-7-2. CASE 7, RUN 1, STEP 2 




Figure IV-7-3. CASE 7, RUN 1, STEP 11 
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. CASE 7, 



Figure IV-7-4 



RUN 1 , STEP 15 




CASE 7, RUN 1 , 



Figure IV-7-5. 



STEP 20 




Figure IV-7-6. CASE 7, RUN 1, STEP 23 
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, RUN 1 , 



Figure IV-7-6P 
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Figure IV-7-7. CASE 7, RUN 2, STEP 1 
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Figure IV-7-8. CASE 7, RUN 2, STEP 2 
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Figure IV-7-9. CASE 7, RUN 2, STEP 9 
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Figure IV-7-10. CASE 7, RUN 2, STEP 12 
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Figure IV-7-11. CASE 7, RUN 2, STEP 14 
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Figure IV-7-12. CASE 7, RUN 2, STEP 27 
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